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The study of production of quark-gluon plasma (QGP) in the laboratory is fascinating for a 
variety of reasons. It is interesting beacuse, this phase did exist for first few microseconds m 
the early universe and perhaps also exists in the stellar matter at high density. Further, a study 
of QGP provides a test of QCD in a new scale (high energy density) in its non-perturbative 
regime. It is also expected to improve our understanding of other fundamental concepts such as 
the nature of confinement and and chiral symmetry. 

From quantum chromodynamics as input dynamics for strong interaction between quarks and 
gluons, we can directly predict this state of matter by large scale computer simulations. The 
results of computer simulations provided a quantitative confirmation of the idea that hadronic 
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matter at sufficiently high density will undergo a transition to a state of matter of deconfined 
quarks and gluons, i.e. quark-gluon plasma On the experimental side, there has been much 
effort m the direction of detecting this phase of matter. In the last decade two experiments at 
SPS of CERN and AGS of BNL were done which succeded in achieving a high energy density, 
although the creation of QGP itself is not confirmed. Future experiments at relativistic heavy- 
lon collidres (RHIC) of BNL and large hadron colliders (LHC) of CERN will hopefully able to 
produce energy density high enough for the production (and equilibration) of QGP 

The most difficult task in these experiments is the detection of quark-gluon plasma. This is 
because, the plasma exists for a very short time in a very small volume and a direct detection is 
not possible The signatures suggested are necessarily indirect and the prominent ones are. (1) 
J /ip suppression, (2) electromagnetic probes such as dilepton and direct photon production, and 
(3) strangeness enhancement. The measurement of these signatures involve both theoretical and 
experimental constraints. The main experimental difficulty lies in singling out the contributions 
arising from the non-qgp sources. The major theoretical difficulty in the calculation of these 
signatures arise from the poorly understood space time evolution of the quark-gluon plasma in 
these ultra relativistic heavy-ion collisions. 

The space time evolution of the quark-gluon plasma evolution at these experiments has been 
generally described by hydrodynamics assuming that the equilibration is almost instantaneously 
reached. In reality, the various stages by which the complete evolution of quark-gluon plasma 
are to be described are, i) pre-equilibrium, ii) equilibrium, where one actually studies the 
thermalised quark-gluon plasma, iii) cooling and iv) hadronisation. It is crucial and very much 
necessary to study the pre-equilibrium stage in the ultra relativistic heavy-ion collisions, for it has 
direct impact on the plasma evolution and also on the measured signatures in the experiments. 
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In this thesis we study the pre-equilibrium evolution of quark-gluon plasma in URHIC by 
employing color flux-tube model. 

In chapter 1 we introduce the topic and provide the necessary background It contains an 
overview of the theoretical and experimental status of the field. A brief description of the 
proposed signatures is presented in this chapter. 

chapter 2 contains an extensive description of the pre-equilibrium evolution of quark-gluon 
plasma within color flux-tube model. In this model the two nuclei that undergo a central 
collision at RHIC and LHC are highly Lorentz contracted as thin plates. These two Lorentz 
contracted nuclei pass through each other and, in the process, acquire a non-zero color charge 
((Q) = 0, ( Q 2 ) 0), by exchanging soft gluons. This leads to creation of a chromoelectric field 
between two receeding nuclei, which polarise the QCD vacuum and produces quarks and gluons 
via Schwinger mechanism (with suitable adaptions and modifications) The quarks and gluons 
so produced collide with each other and get accelerated by the background chromoelectric field. 
In this chapter we set up the transport equations in an extended phase space (which include 
color in addition to coordinate and momenta) to study the production and equilibration of 
quark-gluon plasma m ultra relativistic heavy-ion collisions. 

In chapter 3 we solve the non- Abelian transport equations for qq plasma in SU(2) gauge group. 
In this exercise which is essentially a warm up for the more realistic study which we perform in 
the next chapter, it is shown that the non- Abelian features are crucial and the so called Abelian 
dominance which was studied so far is not justified. We obtain the distribution function of 
quarks and gluons which carry explicit color charges those are abscent in any other study. The 
physical quantities such as local energy density, number density and temperature of the plasma 
are determined for this qq plasma. 
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In chapter 4 we again solve non-Abelian transport equations to include gluons m addition to 
qq m SU(3) group which is the realistic gauge group for strong interaction We made a detailed 
comparison with the SU(2) results as well as with other p-QCD based approaches such as parton 
cascade model. 

A calculation of the dilepton rate is presented in chapter 5 using the distribution of quarks 
and antiquarks within color flux-tube model. It is shown that the pre-equilibrium production of 
dilepton pair is dominant for very small transverse mementum of the dilepton pair, whereas the 
Drell-Yan contribution is dominant over large transverse momentum of the dilepton pair 

In chapter 6, we propose a new source term for quarks and gluons as is appropriate for URHIC. 
We treat the time dependence of the background field dynamically, and also project the amp- 
litudes to a coherent basis in the extended phase space, so as to obtain a legitimate classical 
source term. We show that the ‘source’ is also necessarily a sink both of which exhibit special 
characters in the phase space. We make a detailed comparison with the source term which are 
obtained using rather ad-hoc prescriptions. 

In chapter 7 we have studied J/v> suppressions m a thermally equilibrating quark-gluon plasma 
at RHIC using short-distance QCD. The evolution of the system is taken from parton cascade 
model. The suppression is found to be 100 percent at RHIC energy. A discussion on the validity 
of the models, which describe the evolution of the quark-gluon plasma, is presented from this 
study. 

chapter 8 contains a very brief literature survey of the other models available for the study of 
the pre-equilibrium evolution of quark-gluon plasma. It also contains a very brief description 
of the lattice approaches towards this field. 
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Chapter 1 


Introduction 


In recent times particle physicists and nuclear physicists have evinced a lot of interest in studying 
a new phase of matter, namely, quark-gluon plasma(QGP). It is theoretically established that 
this state of matter consists of unbound quarks and gluons at a temperature about 200 MeV. 
Such a phase did exist in the early stage of the evolution of our universe, ~ 10 -4 second after 
Big Bang. As the universe expanded and cooled down, it went through a phase transition to 
the hadronic matter that we see today. To study this phase and to re-create the early universe 
in the laboratory is a challenging task for both the theorists and the experimentalists. 

Quantum chromodynamics (QCD), which successfully describes the strong interaction between 
quarks and gluons, has predicted the possibility of such a phase. Within QCD, quarks carry 
three types of color charges which interact among themselves via gluons, which are the non- 
abelian gauge bosons. Unlike the Abelian theory where photons do not carry any electric charge, 
gluons in the non- Abelian theory (described by the gauge group SU(3)) carry a charge which 
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is a vector m eight-dimensional space. Within this theory, quarks belong to the fundamental 
representation of SU(3) group and gluons to its adjoint representation This theory explains 
many phenomena observed at the present day experiments. 

However the full range of QCD is not solved till today. This is because, the perturbative 
calculation is not applicable over the large distance scale due to the nature of the strong in- 
teraction between quarks and gluons. This regime needs to be studied within the frame work 
of non-perturbative QCD So far, computer simulations on the lattice[l] provide the only reli- 
able method available for the non-perturbative study. More importantly lattice QCD predicts a 
phase transition between the hadronic state to a quark-gluon phase. It also predicts that at a 
critical temperature T c the dimensionless energy density measure, e/T 4 , increases abruptly from 
a value near that of an ideal pion gas (~ 1) to a value that of an ideal quark-gluon gas. This 
increase in the energy density, Ae, represents some thing like “latent heat of deconfinement”. 
This suggests that for a stable and large system of vanishing net baryon density, there occurs 
a deconfined phase of the quarks and gluons, called quark-gluon plasma. This deconfinement 
phase transition is expected to occur at a temperature of the order T c = 150 A/e V corresponding 
to an energy density of t c — 1-3 GeV//m 3 . Energy densities of this magnitude which did exist 
m the early universe is expected to be obtained in small volumes (~ 100/m 3 ) and for a short 
time (~ 5 — 10/m) in high energy heavy-ion collisions. Abundant multiple scattering among the 
incident nucleons in the heavy-ion collisions and the many secondaries produced during these 
processes are expected to give rise such an energy density. 

As no direct observation exists in these experiments for the initial energy density released, it 
has to be determined from the distribution of the final state hadrons. The formula obtained 
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from the Bjorken model[2] which is generally used to estimated the energy density is 

_ dEr/dy 
nR 2 T 0 


( 1 . 1 ) 


where r 0 = 1/m being the formation time and R = l.lSA 1 ^ 3 . The quantity dEr/dy in the 
above equation is proportional to the final multiple distribution dN/dy which increases as the 
number of binary nucleon-nucleon interaction increases. In order to obtain such an energy 
density through large multiple collisions, it is required that the two nuclei in the collision should 
be large and accelerated to high energy. 


During the last decade such experiments involving heavy-ion collisions have been done at the 
AGS[3] of BNL and the SPS [4] of CERN with Ei a b = 14 and 200 A GeV. Both are fixed target 
experiments and hence produce a low centre of mass energy ( E crn ). Although the results of 
these experiments suggest an extremely rich physics which can not be explained by a straight 
forward extrapolation of nucleon-nucleon collisions, it is still not clear whether (an equilibrated) 
quark-gluon plasma is created in these experiments, in particular at SPS (Pb-Pb collisions at at 
E cm = 17 GeV). The future heavy-ion collider experiments at RHIC (Au-Au collision) of BNL 
and LHC (Pb-Pb collision) of CERN with E cm = 200 A GeVand 5.5 A TeV respectively will 
provide the best opportunity to study such a phase. These values of high energy imply that 
at the moment of maximum nuclear overlap the energy densities are of the order of e, ~ 250 
(7500) GeV/fm 3 at RHIC (LHC), which is much higher than the critical energy density for 
deconfinement caluculated from lattice QCD. 


The most important task on hand is a detection of QGP, assuming that it is produced. As 
the plasma exists for a short time (~ 5 — 10/m) over a small volume (~ 100/m 3 ), it is not 
directly accesible in these experiments. Indirect signatures are proposed to detect this phase 
and the promments among them are : (1) J/ip suppression]/], (2) electromagnetic probes such 
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as dilepton and direct photon production[6, 7], and (3) strangeness enhancement [8]. 

• J/xp suppression 

According to the proposal of Matsui and Satz, the J Jxp production is suppressed in the presence 
of a deconfined QGP medium due to Debye screening. At high temperature the Debye screening 
length which is calculated from lattice QCD is smaller than the J/xp bound state radius. This 
forbids the binding of cc to J/xp. Although suppression of J/xp is seen in heavy-ion collisions, 
still it is not clear that it is due to screening only This is because, J/xp is also supressed 
m p-A collisions and A-B collisions at SPS, where there is no QGP phase. This can be seen 
from Fig-1.1 where the survival probability of J/xp is plotted against various nuclei in different 
experiments. 

The suppression of J/xp m the above cases is due to the presence of a nuclear medium. The 
inelastic collision between cc and nucleon explains the data in p-A collisions and A-B collisions 
at SPS. Although a recent measurement of J/xp by NA50 experiment (see Fig-1.1) is not ex- 
plained by the presence of a nuclear medium, it is still not clear whether this supression can be 
understood in terms of the presence of a deconfined medium. Recently Kharzeev and Satz have 
employed the results of Bhanot and Peskin (1979) to study the interaction of heavy quarkonium 
with gluons. This work provides a dynamical basis to study the J/xp suppression by using 
short-distance QCD. Within this formalism J 'ftp dissociates to D D due to the interaction with 
hard gluons which are present in the deconfined medium. This mechanism is expected to give 
rise a huge suppression of J/xp at RHIC and LHC[9, 10]. Besides this, various authors have 
devloped different proposals for explaining the NASO data without assuming any deconfined 
medium[ll-14]. 
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1/3 1/3 

A + B 

Figure 1.1: The J [rjj survival probability after absorption through nuclear matter. Here A and 
B are the mass number of the colliding objects. The full line is the sutvival probability for 
the proton-nucleus (nucleus-nucleus) systems, calculated with a cross section u = 6.2 mb. The 
marks are the experimental datas, as mentioned in the figure. 


In short, although a supression of J/^ is seen in these existing experiments, no clear conclusive 
evidence for the existence of quark-gluon plasma can be assured. A detailed description of the 
recent theoritical and experimental developements on J/t/i production and its suppression is 
presented in chapter 7. 


• Electromagnetic probes 

Dileptons and single photons have long been proposed as useful probes of the plasma[7]. Once 
produced, they hardly interact with the strong matter and thus carry the details of the circum- 
stances of their production. The major processes for dilepton and photon production in URHIC 
are, (i) hard parton scattering, (ii) electromegnetic decay of hadrons and (iii) production from 
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partons present in QGP (thermal production). The hard parton scatterings produce high p t 
lepton pairs and photons which are calculated from pQCD. The leading order processes are 
quark-antiquark anihillation and quark-gluon compton scattering. This is simillar to hadronic 
collisions. On the other hand, electromagnetic decay of produced hadrons is the main source 
of dilepton and photon production which obscure the signal of interest (thermal production), 
which are produced from QGP. More importantly, the thermal emission of dilepton and direct 
photon is not well understood either theoretically (due to lack of proper description of the space 
time evolution of QGP) or experimentally. These thermal lepton pairs and thermal photons 
are expected to be important signatures of QGP because, their measurement is probably the 
only possibility to have a direct access to the transient dense and high temperature phase of a 
heavy-ion collision. Using thermal parton distributions, the emission rate of thermal dilepton 
and photon is calculated in a simillar way as the high p t direct photon and dilepton calcula- 
tion using the structure functions of nucleons inside nucleus. It is expected that in high energy 
heavy-ion collisions such as at RHIC and LHC, the thermal production will be more than the 
other processes. 

At present, a lot of effort are being made both theoretically and experimentally to understand 
the thermal production of dileptons and photons from the quark-gluon plasma. This is one of 
the hot topic in the search of QGP. A detailed description of this topic is presented in chapter 
5. 

• Strangeness enhancement 

Strangeness enhancement was proposed as a possible signature of the QGP state, almost 17 
years ago[15]. Recall that QGP is formed when the matter density is 10 times greater than 
that of the nuclear matter. The energy corresopnding to the temperatute which is greater than 
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that of the critical temperature is larger than the mass of the ss pair (~ 300 MeV). Therefore, 
the ss will occur at almost the same rate as uu or d d in the QGP at high temperature. In the 
hadron gas the strange particles are produced through associated production during the collision 
of two typically non-strange hadrons The threshold energy to produce this strange particle is 
determined by the mass of the strange hadron pair. The reaction with the lowest threshold is 
p + n— > A 0 + K + + n which requires 671 MeV. On the otherhand. as we have discussed, in 
a QGP the threshold for the strangeness production is only around 300 MeV. Therefore, the 
strangeness production m the QGP is favoured from the production in the hadron gas. The 
strangeness enhancement has been investigated by many experimental groups[16] and recently 
it has been one of the important topics m relativistic heavy-ion collisions. 

The signals discussed above had been initially studied by assuming an equilibrated QGP. 
However, the measurement of these signals in the actual experiment, involve the contributions 
from all the stages of evolution of the system. In order to make a convincing case for observation 
of QGP, the observation of several signals must be supported with a reasonable description of 
the interactions of the final state particles. The formation and study of quark-gluon plasma 
in future experiments, will crucially depend upon a detailed microscopic understanding of the 
plasma. Questions relating to its formation, evolution and hadronisation deserve a much better 
study. 

Of particular importance to us is the production and the evolution of quark-guon plasma after 
two nuclei collide with each other. It is expected that m the center-of-mass frame, two highly 
Lorentz contracted nuclei pass through each other and deposit their kinetic energy gradually 
behind them[17]. A qualitatively accepted space-time evolution of these reactions is schematic- 
ally illustrated in Fig-1.2. The various stages by which the complete evolution of quark-gluon 
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Figure 1.2: A typically accepted space time evolution picture of the nucleus-nucleus collisions 
in ultra relativistic heavy-ion collisions 


plasma is described in URHIC are, 1) pre-equilibrium, 11) equilibrium, where one actually stud- 
ies the thermahsed quark-gluon plasma, in) cooling and iv) hadronisation. In the very early 
stage of this collision, the partons are being produced and they collide with each other m the 
pre-equilibrium stage to form the equilibrated quark-gluon plasma m the stage (ii). It is expec- 
ted that most of the entropy and transverse energy are produced in this satge. The dynamics 
of partons during the first few fermi determines the initial (and boundary) conditions for the 
hydrodynamic expansion of the plasma phase. The validity of the calculations involving hy- 
drodynamic expansion, assuming very fast local thermodynamic equilibrium after the partons 
are produced, is to be checked and a proper study of pre-equilibrium dynamics is needed. In 
particular, we want to know how the partons evolve, under what conditions they thermalise, 
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the time scale of equilibration and the typical energy density soon after the thermalisation 
Without these inputs there are always uncertainties in choosing the equilibration time and the 
initial temperature which crucially determines the predictions of all the signatures. To study 
these questions systematically in a realistic model, a detailed study of the pre-equilibrium stage 
during the very early stage of the high energy nuclear collisions is necessary. 

As mentioned, the production of an equilibrated quark-gluon plasma in stage (li) crucially 
depends on the pre-equilibrium evolution, and hence all the signatures in the experiments. 
Indeed it has been pointed out [9. 10] that J/ip suppression which is generally attributed to the 
equilibrium stage can also arise from the pre-equilibrium stage, an observation which is true for 
strangeness enhancement as well. The contribution of the pre-equilibrium stage to the dilepton 
production and the direct photon production would be even more pronounced, and a careful 
study of the above signals ought to shed light on this stage of the plasma. 

In this thesis we discuss the production and the equilibration of quark-gluon plasma(QGP) m 
URHIC by employing the Flux-tube Model[18, 19] which is a generalization of the familiar Lund 
string model[20] widely used for e + e~ and p — p collisions. A chapter wise summary is given in 
the synopsis (see page v). 



Chapter 2 


Pre-equilibrium evolution of quark-gluon 
plasma 


2.1 Introduction 


The space-time evolution of quark-gluon plasma created in ultra relativistic heavy-ion collisions 
is governed by hydrodynamics, once it reaches equilibrium The relativistic hydrodynamic 
equations are: 


and 


dj^{r) = 0 

(2.1) 

c^n M (r) = 0, 

(2.2) 
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where T M "(r) is the local energy-momentum tensor and n^(r) is related to the local number 
density n(r) by n(r) = n M (r)u M (r) with u M (r) being the flow velocity. In the absence of any 
dissipation the energy-momentum tensor is written in terms of energy density e(r) and the 
pressure p(r) as 

T M1/ (r) = -g^pir) + (e(r) + p(r))u M (r)u 1/ (r). (2.3) 


where e(r) = T M1/ u M (r)u l/ (r). If one further neglects the transverse expansion in the early stage 
and considers only the longitudinal expansion, we then arrive at the equation[2]: 


dc t + p 

dr t 


= 0 , 


(2.4) 


first proposed by Bjorken to describe space-time evolution of QGP In the above equation, 
r = J ( t 2 — z 2 ) is the proper time. 


As these hydrodynamical equations are applicable only in equilibrium, which is attained at 
certain time after the two nuclei have collided with each other it is desirable to determine 
the equilibration time by some microscopic approach. After all, the initial conditions that 
govern the hydrodynamical evolution (Eqn 2.4) can only be obtained from a careful study of 
the pre-equilibrium regeime in heavy-ion collisions For an ideal quark-gluon plasma of massless 
partons, the equation of state is given by [17] 


p = a(T)T 4 


and 

e ~ 3 a{T)T 4 (2.6) 

s = ilZ ~ 4a(T)T 3 (2.7) 

and 


n = 6(r)T 3 . 


( 2 . 8 ) 
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Here s is the entropy density, T is the temperature, n is the number density of the plasma and 
the quantities are given by 


7T 2 21 

a ( T ) = qq( 16 + y n f) 

(2.9) 

b(T)= ^(16 + 9"/) 

'77* “ 

(2.10) 


with 7 if being the number of flavours. It follows from the above equation of state together with 
Bjorken’s scaling picture that the temperature of the plasma scales like T(t) cc r~C 3 . energy 
density like e(r) a r~ 4 / 3 and number density like n(r) cc t~ x . 

The major difficult task is the detailed study of the poorly understood pre-equilibrium stage of 
the collisions Crudely speaking, this pre-equilibrium stage is the duration between the moment 
when the two nuclei just pass through each other to the actual establishment of an equilibrated 
quark-gluon plasma. While most of the earlier work refers to the equilibrium stage studied by 
employing the hydrodynamic evolution, it is necessary to study the pre-equilibrium stage for 
the above mentioned reasons. Indeed, without a proper study of the pre-equilibrium stage there 
are always uncertainties in the predictions of all the signatures in the ultra relativistic heavy- 
lon collisions (at RHIC and LHC). Any pre-equilibrium study also presupposes a production 
mechanism for QGP, and here, we study the pre-equilibrium evolution of quark-gluon plasma 
in ultra relativistic heavwon collisions (URHIC), within color flux-tube model. The color flux- 
tube model is a generalization of the familiar Lund string model widely used for e + e~ and p — p 
collisions[20]. A brief description of the model is given below. 
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2.2 Color Flux- tube model 

In this model the two nuclei that undergo a central collision at RHIC and LHC are highly lorentz 
contracted as thin plates[21, 22] (The two nuclei in collision at RHIC and LHC travel almost 
at a speed that of light). These two lorentz contracted nuclei pass through each other and, in 
the process, acquire a non-zero color charge (( Q ) = 0, (Q 2 ) ^ 0), by exchanging soft gluons. 
So one may figuratively call such nuclei after collision as color capacitor plates between which a 
strong chromoelectric field is created. The nuclei which act as color capacitor plates produce a 
chromo-electric field between them [21, 22] The strength of the field which naturally depends on 
the strength of the color charge residing on the plates cannot be fixed from first principles. We 
can only fix that phenomenologically, say by identifying the field energy with the energy in the 
central region as estimated by Bjorken[2] This strong electric field creates qq and gluon pairs 
via the Schwinger mechanism which enforces the instability of the vacuum in the presence of 
an external field. The partons so produced, collide with each other and also get accelerated by 
the background field. The mutual collisions drive the system towards equilibrium with suitable 
modulation from the background acceleration Within this model a set of transport equations 
are solved to study production and equilibration of quark-gluon plasma in URHIC. A brief 
description of works done on this model is given below. 

In the approach of Baym[l8], the partons are assumed to be created soon after the two nuclei 
have collided with each other, and are further assumed to equlibrate almost instantaneously 
thereafter (with a very small collision time). This quark-gluon gas then proceeds to expand hy- 
drodynamically. Subsequent work of Kajantie and Matsui[19], incorporated a dynamic particle 
production by introducing a source term via the Schwinger mechanism[23] in the Boltzmann 
equation, and they also studied the equilbration in the relaxation time approximation. On the 
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other hand, Bialas et a/[24] generalized the Baym analysis by introducing the effect of the back- 
ground field on the otherwise hydrodynamic flow. A self consistent study of the system was 
carried out by Banerjee et al[ 25] who combined the effect of both the background field and the 
collisions on the quark-antiquark system. Finally, by employing the same analysis as of Banerjee 
et al, Asakawa and Matsui[28] have studied the rate of dilepton production by qq annnihilation. 
as a function of (proper) time All the above analyses ignore the gluon component. 

However some of the crucial features which have not been considered m these earlier studies 
remains to be justified. First of all m the analyses of Kajantie et a/[19] and Banerjee et a/[25], the 
source term acquires a time dependence by virtue of the time dependence that the electric field 
suffers because of the particle production The production of particles is then strictly go\erned 
by a time dependent field, in which case a perturbative mechanism will take over from the non- 
perturbative Schwinger mechanism. Note that the latter holds only for a uniform constant field. 
All the previous studies employing the flux tube model employ the non-perturbative expression. 
2) The source term is derived for QED processes and needs to be redone for QCD processes at 
hand It is apparent that the QCD effect will show up most manifestly in the gluon production; 
there is no corresponding counterpart in QED. In other words it is necessary to obtain a source 
term for gluons. 3) It is well known from earlier studies of Yang-Mills equations that there 
are inequivalent gauge field configurations which yield the same field tensor[26], and hence the 
same energy momentum tensor. This feature which is again absent m the Maxwell case has to 
be taken into account m this study. Recall that the lattice studies at any given energy density 
implicitly sum over the contributions coming from all such inequivalent configurations^?]. 4) 
Finally, the inherent non-linearity and the local gauge invariance which characterize the Yang- 
Mills equations lead to an added complexity; the equations do not possess a unique solution, 
even after gauge fixing, and even if one considers only static configurations. Consequently, one 
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needs to scan, with some suitable weightage, all the solutions for the given initial conditions, 
and sum over all such configurations It is not clear if any procedure is known that allows an 
implementation of this requirement 

We have attempted here to fill this gap partially within color flux-tube model by solving transport 
equations m an extended phase space 


2.3 Transport equation in an extended phase space 


In the study of production and equilibration of quark-gluon plasma in ultra relativistic heavy-ion 
collisions one needs to know the distribution of quarks and gluons, f{x,p ), in the usual phase 
space. The determination of all physical quantities in the experiments crucially depend on this 
distribution functions of quarks and gluons. At equilibrium the above distribution functions 
are the standard Fermi-Dirac distribution function for the quarks and anti-quarks and Bose- 
Einstein distribution function for gluons. In the non-equilibrium stage the time devlopement of 
the phase-space distributions f a (x,p ) for the partons of species a = qj,qj, g (f being the quark 
flavors) is described by the relativistic transport equation (Boltzmann equation): 


P nd 11 + fa{x,p) = Ca{x,p) + Sa{x,p) 


( 2 . 11 ) 


with some known initial condition of f a (x,p) at t — t 0 , when the two nuclei begin to overlap at 
RHIC and- LHC. The first term in the left hand side of the above equation corresponds to the 
usual convective flow and the second term is the Lorentz force term. 5 and C on the right hand 
side correspond to the source term for parton production and collision term among the partons 
(described in chapter 3). 
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The above equation is written in the usual six dimensional phase space and does not carry 
any information on the time evolution of the color charge, which is a dynamical quantity for 
non- Abelian plasma. Thus, as such, these transport equations which are valid for an abelian 
plasma are not applicable to quark-gluon plasma. It is our purpose to emphasize and study a 
crucial feature of the quark-gluon phase. viz, the inherently non- Abelian nature of the plasma. 
(Note that all the earlier macroscopic approaches[18. 19, 24, 25, 28] are in the so called Abelian 
approximation). Put simply, the non- Abelian features are completely ignored m the space inme 
evolution, and the Boltzmann equation describes the dynamics which is essentially that of a 
Maxwell fluid, be it in the source term or in the background field term. The approach of Geiger 
and I<apusta[17] does incorporate the non-Abelian features in their basic Feynman diagrams, 
but is neverthless incomplete m that the distribution function in its final form does not carry 
any color degee of freedom. 

We will now set up the relevant equations with the color degree of freedom, following several 
works of Heinz[29]. In non-abelian theory the color charge is a continuously varying function of 
time. The precession of the color charge(Q a ) obeys Wong’s equation[30]- 


rlD a 

— = f ahc u^Q b A ctl 
dr 


( 2 . 12 ) 


which supplements the Lorentz force equation 


djA 

d7 


Q a F aiiU u u . 


(2.13) 


Here A a M is the gauge potential, and f ab ~ is the structure constant of the gauge group. In older 
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to include the color charge m the phase-space, we consider an extended one particle phase space 
of dimension d — 6 + (/V 2 — 1) in SU(N) The extended phase space is taken to be the direct 
sum R 6 © G, where G is the( compact) space corresopndmg to the given gauge group In short, 
in addition to the usual 6 dimensional phase space of coordinates and momenta we now have 
another eight coordinates corresponding to the eight color charges in SU(3). In this extended 
phase space a typical transport equation reads as[31] - 


+ Q a F:yd£ + r bc Q a Ay-d c Q 


f(x,p,Q) = C(x,p,Q) + S(x,p,Q) 


(2.14) 


Here f(x,p,Q ) is the single particle distribution function in the extended phase space The 
above equation is gauge and Lorentz m\ariant[31]. The first term m the left hand side of 
equation (2.14) corresponds to the usual convective flow, the second term is the non- Abelian 
version of the Lorentz force term and the last term corresponds to the precession of the charge 
as described by Wong's equation. Note that, in general we have to write separate equations for 
quarks, antiquarks and gluons since they belong to different representations of the gauge group. 
For anti-quarks the distribution function /(x,p, Q ) obeys a simillar equation, with Q a replaced 
by — Q a ( i.e . the second term in the above equation changes sign). These equations are closed 
with the Yang-Mills equation, given by 

[D^F liU ) a {x) = /(x) = g J p l/ Q a [f q (x,p,Q) - f q (x,p,Q) + f g {x.p,Q)]dPdQ (2.15) 
The integration measures in the above equation are 

£ (P=2e(p„)i( P 2 -m 2 )4Y = -^|- (2.16) 

(2tt ) j {ZttFpo 


dQ = 8{Q a Q a - q 2 )5(d abc Q a Q b Q c - f)d 8 Q. 


(2.17) 
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Using these distribution functions we obtain the local number density and local energy density 
of the quark-gluon plasma by the equations 

n (r) = y dPdQip^u^if^x.p.Q) + f q {x,p,Q) + f g {x,p,Q)) (2.18) 

e(r) = J dPdQ(p»u^) 2 {f q (x,p,Q) + f q {x,p,Q) + f a {x,p,Q)) (2.19) 

In the next chapter we shall consider the simplest of the non- Abelian group viz SU(2) to highlight 
the dynamical role of the color degree of freedom. The realistic guage group SU(3) applicable for 
QGP will be studied m detail in chapter 4 Both the chapters contain rather ad hoc adaptation 
of the Schwinger formola for qq and gg production that will be reconsidered in chapter six. A 
full study of the transport equation with the new source term is beyond the scope of this thesis 
That will be reported elsewhere. 



Chapter 3 


Solution of transport equations and 
evolution of qq plasma 

3.1 Transport equations for q q plasma 

In this chapter we discuss the production and the evolution of quark-antiquark plasma in URHIC 
by employing the Flux-tube Model[l8, 19], in particular to obtain a time devlopment of the quark 
and antiquark distribution function. As mentioned earlier, m this model the two nuclei that 
undergo a central collision at ultra high energies are lorentz contracted as thin plates [21, 22]. 
These two lorentz contracted nuclei pass through each other and, in the process, acquire a non- 
zero color charge ((Q) = 0, ( Q 2 ) ^ 0), by exchanging soft gluons. So one may figuratively call 
such nuclei after collision as color capacitor plates between which a strong chromoelectric field 
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is created. This field creates qq and gluon pairs via the Schwinger mechanism which enforces the 
instability of the vacuum m the presence of an external electric field The partons so produced 
will collide with each other and also get accelerated by the parent background field. The mutual 
collisions drive the system towards equilibrium with suitable modulation from the background 
acceleration. In all the earlier studies (e g [19, 25]) a set of transport equations are solved within 
this model in the abelian approximation, i.e. the model was valid for an electromegnetic plasma. 
We now attempt to study the evolution of a non-Abelian plasma by solving transport equations 
which are written in an extended phase space of coordinate, momenta and color. 

As mentioned eariler, for simplicity and also to get a feeling for how crucial is the non-abelian 
effect in the study of quark-gluon plasma can be, we shall study the gauge group SU(2). and 
consider only a qg plasma. We will solve the transport equations in SU(3) in the next chapter 
with gluons taken into consideration. For SU(2) charges that we are interested in. the non- 
Abelian extension of Eqn.(2.14) reads 


+ Q a F;yd; + e abc Q a A b yd c Q 


f(x,p,Q) = C{x,p,Q ) + S(x.p,Q). 


( 3.1) 


Note that we do not have to write a separate equation for anti-quarks since they belong to the 
same color representation as quarks. The quarks considered m this chapter does not possess any 
other quantum number such as the electric charge. In reality the situation is of course different. 
Indeed, the antipodal points on the sphere (corresponding to the color part of the phase-space 
with a fixed charge) represent the particle and anti-particle. Wong’s equation (Eqn 2.12) may 
be recast in the more convenient form 
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dO 

dr 


= x ^ 


(3.2) 


where the arrows now denote the direction m the color space. In the same notation, the transport 
equation obtains the form 


f(x,p, Q) = C(x,p,Q) + S(x,p, Q). (3.3) 

For the model at hand, let us take the "plates’ to be moving along the z-direction. Also we 
restrict to only that field configuration for which there exists a gauge choice such that the 
gauge potentials commute with each other every where. Keeping this m mind, we make the 
gauge choice where only the components A M “ = (A 03 , Al 33 ) are non-vanishing and the other 
components are zero. With this choice, the resulting electric field points in the ‘3’ direction in 
the color space. We next require a boost invariant description[2] of the distribution functions 
as well as the other physical quantities which may be determined thereof, even as the system is 
evolving. We now proceed to impose the Lorentz gauge condition, which has to be done such 
that the chromoelectric field depends on the boost-invariant(along the axis of collision) quantity 
t — (t 2 — z 2 ) 1 ! 2 Observing that the system is effectively (1 + 1) dimensional, m the t-z plane, 
we write 


+ Q ■ F^d; + p^Q X A u - 


_d_ 

dQ 


A* = e^d„G{r) (3.4) 

where the indices p, u take values 0, 3 and G is a Lorentz scalar function depending only on r. 
The above choice automatically implements the Lorentz gauge condition = 0. Clearly, the 
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chromoelectric field E is dependent only on r, and is given by 

G(r) (3.5) 

(It will be seen below that all the derivation in this chapter is obtained in terms of E(r), not 
G(r). However this relation (Eqn. 3.5) is important (see chapter 6) when we use new source 
term (from r dependent electric field) for particle production.) 

The Wong equation (Eqn 3.2) guarantees the conservation of the magnitude of the vector charge 
Q, which may now be held fixed. Being the analogue of the Larmor equation for a charged 
particle in an external magnetic field, it also conserves the component of the charge that is 
parallel to the background chromoelectric field. It is, therefore, convenient to resolve the SU(2) 
charge in the polar coordinates. Writing 


E(r) = 


dr . . d 
+ 


Q i = QsmOcoscf ) ; Q 2 = QsinOsiruf)-, Q 3 = QcosO , 


(3-6) 


it is straight forward to verify that 


Q X dQ) 3 d<p 


(3.7) 


The last equation will lead to considerable simpification in solving the transport equation. 


Let us now consider the collision term. The color plate model can be looked upon as an effective 
version of a more complete microscopic theory describing the interactions of quarks and gluons. 
Hence, we employ the simple relaxation time hypothesis, with a phenomenological parameter 
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r c , the relaxation time, to be obtained from microscopic computations. It then follows that 


C = 


~PX(/ Z Jjll 

r c 


(3 8) 


where f eq is the eqilibrium distribution function, with local (space time dependent) values of the 
thermodynamic quantities. While it has been customary to take f eq to be one of an ideal gas, 
with some support from the parton cascade model (pcm) analysis of Geiger and Kapusta[17], 
recent lattice studies suggest that the quark-gluon phase is possibly a non-perturbative phase, 
being ideal only at T ~ 5 T c . Since we are not dealing with a true system here, we shall 
conveniently take f eq to be that of an ideal Fermi gas with a local temperature, evolving as 
a function of the proper time r Apart from space time dependence this distribution function 
carries color in its definition [32]. We thus have 


/e 9 = 


2 

exp((p^ - A fl )u fl /T(r)) + 1 ’ 


where is the flow velocity, 


(3.9) 


u M = (coshrj, 0, 0, sinhrj). (3.10) 

which is written m terms of the (space-time) rapidity tanhr] = z/t. 

It can be verified that with the above choice of the potential 


AX = 0, 


(3.11) 
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and the above equation boils down to the familiar form 

2 

^ exp (p*)uJT(r)) + 1 

Now we demand the boost invariance following Bjorken’s picture according to which the lon- 
gitudinal boosts are the symmetry operations on the single particle distribution. The boost 
invariant parameters on which / can depend are, apart from the charge coordinates, 

r = (t 2 - z 2 ) 1/2 , £ = (77 - y), p t = (Po ~ Pl) 1/2 (3-13) 

where y = tanh -1 (p;/po) is the momentum rapidity. It is convenient to write separate equations 
for quarks and anti-quarks, a la the Abelian case. If we therefore identify the quark states 
with the points on the upper hemisphere of the color sphere, and the antiquarks with the lower 
hemisphere, by a trivial relabeling, we may write two equations, 


(3.12) 


d_ 

dr 



+ 


g cos 9E(t) 
Pt cosh £ 


d 

dt 


+ g—G(r) tanh £ 

dr 


d 

d4> 


/(r,£,pi,0.6) 


+ L _ hi + 

r c r c p t cosh £ 


(3.14) 


d ( tanh £ 

dr \ r 


g cos 9 E(t)\ 

p t cosh £ J 


d 

3£ 


d 

+ g—G{r ) tanh£ 

ar 


d 




+ L _ hi + 

r c T c p t cosh £ 


(3.15) 


with the first of them for the quarks and the next for the anti-quarks. The angle variable 9 
varies from 0 to tt/2 in both the equations. 
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Finally, S is the (non perturbative) Schwinger’s expression for pair production and is given by 


S(r,£,p t .6>) = 


gE cos 8 
8tt 3 


In 


— exp 




gE cos 8 


(^) 1/2 exp (— a£ 2 ) 


(3.16) 


where we have inserted the Gaussian dependence on £ by hand. Note that the non-Abehan 
nature features dominantly via the cosd term This occurence plays an important role m the 
evolution of the system. 


After having set up the relevant equations (not all yet since energy-momentum conservation is 
to be imposed), we observe that the above differential equation possesses the formal solution 


f(T,Z,Pt8,4>) = / dr exp ( 
Jo 



s (T'^Pt.fl) | Ur^Ptl 

Pt cosh E T c 


where £(r/) is given by 


(3.17) 


= sinh 1 


— sinh £ + 


g cos 6 

Pt-r' 



c1t"E(t ") 


(3.18) 


It is thus clear that / does not depend on o, so the term contributes nothing to the transport 
equation. However the 9 dependence is still involved through out the formalism which shows 
the non-Abehan effects. The corresponding distribution / for antiquark can be written by just 
changing g to —g m / as can be checked out easily 


The transport equation has to be supplemented with the conservation constraint 
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d.TZt + d.TZr = 0 , 


(3 19) 


since it is the field energy that is being pumped in order to produce the qq pairs More explicitly, 


ajz, = -fjr = ( 3 . 20 ) 

with 

TZ, = Jp“p“(f + f)dVdn Q (3.21) 

where the measures cLVLq = sin 0d9d(f> , dT = an d 3a = f P^Qaif ~ f)dVdflQ. 

Here we have taken the value of the degeneracy factor 7 = 2 corresponding to two flavours. 
Since energy and momentum are conserved in each collision, the moment of the sum of the 
collision term vanishes. 


jf(C)dTd9. Q = 0 (3.22) 

Now taking the first moments of the Boltzmann equation and integrating over the color degrees 
of freedom for / and / and making use of the conservation of energy and momentum, we obtain 


djr + gE{r) dTdV. Q p' 


8<J - /) 

of. 


+ 2 / dVdilgS = 0 


(3.23) 


where 
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Tj u = diag(£ 2 /2, £ 2 / 2, £ 2 /2, -£ 2 /2) (3 24) 

is the energy momentum tensor for the field We may solve for the electric field by employing 
the same procedure as m[25], and by employing the symmetry /(r, f.px,#) = /(t, — f,px, 0). 
We thus obtain the equation governing the decay of the source field to be 


dEjr) 

dr 


Jo dptp2t Jo d Z smh -tJ 0 d(9sin6>[/ - /] + ya|£(r)| 3/2 = 0 


(3.25) 


where a = a£(5/2) exp(0.25/a), a = c(g/ and c = jj. Finally, q (5/2 ) = 1.342 is the 
Reimann zeta function. 


Equations (3.25) and (3 17) are as yet underdetermined since the local temperature T(r) is free. 
In order to fix the form of T(r), we appeal to the relaxation time approach that we are employing, 
and assume that the the particle energy density differs negligibly from the equilibrium energy 
density. We may then relate, by an ansatz, the proper energy density which is defined by 


e(r) = j drdn Q (p“u r ) 2 (f + f) (3.26) 

to the temperature by its equilibrium value, whence, 

TW = ((7^3) £( r))‘ /4 ' (327) 

It may be mentioned that the weaker condition of energy-momentum conservation that we have 
employed here , m fact, satisfies the Yang-Mills equations D P .F^ ,L ' = j' J as well. 
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Finally, we pause briefly to discuss the effects of hard thermal loops that have been emphasized 
recently [33]. They have been derived from the classical transport equation as well by Kelly 
et al{ 34, 35] The latter derivation, which is important for us here, is based on an analysis of 
the Vlasov equation which properly describes the expansion of an already equilibrated gas, m 
a back-ground field but with with no source or sink. It is. therefore, necessary to determine 
whether the more complete transport equation such as the one that we have here will have any 
effect on the results of Ref. [34, 35] Conversely, it is also necessary to determine how the hard 
thermal loops will affect the purely classical results that will be obtained here. Let us recall that 
the derivation of Kelly et al[34, 35] consists of a systematic expansion of the distribution function 
as well as the background field term in powers of the coupling constant. The zeroth order term 
is merely the Fermi-Dirac term for the quarks. Since we are solving here the approach to the 
equlibrium, there is necessarily a dependence of the distribution functions and other physical 
quantities on g. It remains to disentangle the contribution coming from the hard thermal loops 
Indeed, observe that the two extra terms that we have at hand here are of higher order in g. 
First of all, the source term is non-analytic and does not even admit an expansion in powers 
of g The other collision term is easily seen to be of order g 4 or higher. We thus conclude that 
the conclusions of Ref. [34, 35] remain unaffected, and we may take over those results in toto , 
supplementing our classical results 


3.2 Computational Procedure 


We have adopted here a double-self consistent method to determine (/, e, ...) following the work 
of Banerjee et al[25] in the Abelian case. The procedure follows the scheme {T(r) tria i, E(r) trta i} — »■ 
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{/, / E(t)} ->{/,/} — >■ T(r) -»... by repeated use of equations (3.17), (3.25), (3.17), (3.27). 
The iteration terminates as soon as a convergence is established in the solutions for E(t).T(t) 
All the desired quantities are thereby consistently determined. The same procedure will be 
employed subsequently. 


3.3 Results and discussions 

Before we proceed to present and discuss the results, a few comments about the choice of 
the value of the parameters in the computation We put g — 4 throughout our calculations. 
Since lattice computation results predict a phase transition from the baryonic phase to the 
quark-gluon phase at densities ~ 5 — 10 GeVj fm 3 , it has also been customary to take an initial 
energy density in the same range m the Flux tube model. However, greater care needs to be 
taken before the initial energy densities in URHIC are chosen. Indeed, a fairly reliable estimate 
of the time required for the partons to be produced in central collisions after the two nuclei have 
suffered the maximum overlap is ~ 0.05 —0.1 /m; a simple dimensional analysis leads to a value 
of the initial energy density to be [1/2)Eq s; 500 — 1000 GtVj fm 3 , the precise value depending 
on the the magnitude of g in \fgE — A Since in the color plate model, we set the zero of 
time not at maximum overlap, but at the instant when the plates have completely crossed each 
other, we shall take the initial energy density e = 300 GeV/ fm 3 . We shall also be guided by the 
results of Geiger and Kapusta[17j m our choice of the values of r c , and take r c = 0.2/m. to be a 
realistic value[17]. The scale for the hydrodynamic limit will be set by the choice for r c as well 
as the formation time r 0 . A typical value for the analysis at hand is 0.001 fm, which we employ 
here. Finally, we study the other extreme case, the collisionless limit, by the choice r c = 5 fm. 
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It is instructive to compare how the realistic regime behaves in relation to the results that we 
obtain for the two extreme limits. 

We have studied the decay of the chromoelectric field, the evolution of the particle energy and 
number densities, the evolution of the distribution function and its approach to the equhbrium 
state We have compared them with the corresponding Abelian results. In addition, we have 
also evaluated a qumtessentially non-Abehan quantity, viz, the expectation value of the angle 
that the quark charge makes with the direction of the chromoelectric field. In the current model 
the chromoelectric field not only decays to produce the qq pairs, but is itself built up by the 
receding nuclei which act as color plates. We have, therefore, also studied the particle energy 
per unit transverse area as a function of ordinary time to highlight this feature. 

3.3.1 Discussion of the results 

We shall present the results of our analysis m the three regimes corresponding to the hydro- 
dynamic, the realistic, and the collisionless cases. Comparison between the Abelian and the 
non- Abelian systems will be taken up subsequently. 

Consider the hydrodynamic limit first. The physical quantity of utmost importance is the 
chromoelectric field, whose dynamics is most readily determined m the Flux tube model at 
hand Indeed, apart from determining the production and the acceleration of the partons m the 
Flux tube model at hand, it has an additional role which has been emphasized by Svetitsky[ 36 ]: 
the charms which are produced in the pre-equilibrium stage would not only interact with the gas 
of light quarks and gluons, before either forming a J/ip or an open charm mesonic state, but will 
also be influenced by the background field. Thus the study of the evolution of the mean electric 
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field m the pre-equilibrium stage acquires an added importance. It is a merit of the Flux tube 
model that we can readily determine the evolution of the background electric field. Note that in 
contrast, no such information can be extracted in the more microscopic models such as parton 
cascade model (PCM)[17] and heavy-ion jet interaction generator (HIJING)[37j. In Fig 3 1 the 
decay of the field is shown Recall that r c = .001/m As can be seen from Fig 3 1. the electric 
field has hardly decayed at all, with a percentage decay less than 2 %, even at r = 1 5/m. 
The corresponding particle energy density, shown m Fig. 3.2 is also negligibly small, with the 



Figure 3.1: Decay of the chromoelectric field as a function of proper time (in units of fermi), for 

t c = . 001/m The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

ratio with respect to the initial energy density being ~ 10~°. The Abelian situation is an 
order of magnitude better, but is still hopelessly small. And indeed, the number density is also 
very small, and as shown in Fig. 3.3 gets saturated at ~ .05/ fm 3 and the temperature also 
stabilizes to a value ~ 40 MeV, which is much less than the temperatures required. Clearly, it 
is very unlikely that the plasma does not go through a prequilibrium phase. 


The behaviour of the system for r c = 0.2 and 5/m is in sharp contrast to the case discussed 
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t (fm) 

Figure 3.2: The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = 001 fm. The solid line refers to the non-Abelian 
case, and the broken line to the Abelian case. 

above There is a significant decay of the electric field, which is shown in Figs. 3.4 and 3.5 
for the two respective values of r c . The electric field decays by about 15% of the original 
value, at r = 1 5/m. The corresponding energy densities are also not very different, and 
yield ~ 10 — 15% of the initial field energy density (see Figs. 3.6 and 3.7). The corresponding 
temperatures are ~ 300A7eVh which is quite realistic. Note that the parton cascade model[17] 
predicts a value about 300 MeV at 2.4/m. The real difference between the collisionless case 
and the ’’realistic” case is m the number density. Indeed as Figs. 3.8 and 3.9 show, the plasma 
produced for r c — 0.2/m has a number density ~ 15 — 20/ fm 3 , the coriesponding number for 
t c = 5/m is three to four times smaller. In other words, the system is more dense in the former 
case than the latter. Clearly, this distinction should show up in signatures such as dilepton 
production which are sensitive to both the number as well as the energy density. 

It thus follows that the flux tube model, with the incorpration of the color degrees of free- 
dom, definitely rules out instantaneous equilibration as envisaged in the approach of Baym[18]; 
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Figure 3.3: The particle number density as a function of proper time (in units of fermi), for 

r c = .001/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

further, it also distinguishes the collisionless limit from realistic values of equilibration time. 

Finally, it should be emphasized that the flux tube model as employed here not only pumps 
in the field energy to particle production, but also contributes to the field energy by virtue 
of the recession of the color charged nuclei from each other. In fact, a crude estimate of the 
time required to convert all the plate energy to the field energy turns out to be ~ 5/m for 
200(?eVy nucleon In any case, it is therefore misleading to interpret the ratios we have shown in 
Figs. 3. 2, 3. 6 and 3.7 as the fraction of the total field energy that has gone into the particles. To 
emphasize this, we have evaluated the dependence of the field and particle energies per unit area 
by integrating over the contribution along the longitudinal direction. The results are shown m 
Fig. 3.10, from which it is clear that there is a lot more energy m the field than in the particles. 
For the same reason, it is also misleading to interpret the field energy density at r — 0 as the 
counterpart of the energy densities employed m lattice computations to study the transition 
from the hadronc to the QGP phase. 
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Figure 3.4: Decay of the chromoelectric field as a function of proper time (in units of fermi), 

for t c = 2/m. The solid line refers to the non-Abehan case, and the broken line to the Abelian 
case. 

3.3.2 Comparison with Abelian results 

It is clear from Figs. 3. 1-3.9 that the Abelian and the non- Abelian results bear little resemblence, 
belying the expectation that the ’’Abelian dominance” holds m this case The difference is most 
dramatically highlighted at r c = 0 2/m. Whereas almost all the initial field has decayed in the 
Abelian case, only 30 7c has done so in the non-Abehan case. Accordingly, The particle energy 
density is larger by a factor of ~ 5 , the number density by a factor of ~ 4, and yields an 
abnormally large value of~ 800 A/e V" for the temperature. Earlier Abelian calculations[19, 25, 
28] yielded reasonable values for the temperature because of unrealistic values for the initial 
field energy It is a universal feature that the non-Abelian plasma is rarer and cooler than its 
Abelian counterpart. In the hydrodynamic limit, the Abelian ananlysis yields a temperature 
~ 200 MeV, in contrast to ~ AQMeV m the non-Abelian case. Thus the incorporation of the 
color degree of freedom rules out instantaneous hydrodynamic evolution. Of course, a colorless 
plasma does not have reasonable temperature for any other r c . 
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Figure 3.5: Decay of the chromoelectric field as a function of proper time (in units of fermi), 

for r c = 5/m. The solid line refers to the non- Abelian case, and the broken line to the Abelian 
case. 

Yet another interesting aspect that emerges from our studies is in the close interplay between 
the value of r c and the color degree of freedom While the field decays faster in the non- Abelian 
case in the hydrodynamic limit, the trend reverses at r c = 0.2/m and gets restored in the 
collisionless limit. In contrast, the particle energy density is smaller for a colored plasma in 
hydrodynamic limit. Further, it continues to be so. at r c = 0.2/m but becomes larger than the 
Abelian case m the collisionless limit. These non-trivial manifestations of the color charge may 
be expected to have an important bearing on the other bulk properties of the plasma. 

In order to gain some insight into the above features, we have also evaluated the expectation 
value of 0 =< cos 2 9 > which yields the rms value that the particle charge makes with the field 
direction ip the color space. Note that the above quantity is gauge invariant, and hence physical. 
Fig. 3.11 shows 0(rj for r c = 0.2/m. It may be seen that the value saturates around 0.25, 
corresponding to 9 ~ r/3. The corresponding value in the non- Abelian case is strictly zero. 
The effect of the background field is thus reduced, leading to dominance of the collision 
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Figure 3 6: The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (m units of fermi), for r c = 2/m. The solid line refers to the non-Abelian case, 
and the broken line to the Abelian case. 



x (fra) 

Figure 3.7: The particle energy density scaled w r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = 5/m. The solid line refers to the non-Abelian case, 
and the broken line to the Abelian case. 
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Figure 3.8: The particle number density as a function of proper time (in units of fermi), for 

r c = 2/m. The solid line refers to the non- Abelian case, and the broken line to the Abelian 
case. 



t(fm) 

Figure 3.9: The particle number density as a function of proper time (in units of fermi), for 
r c = 5/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian case. 
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term in the expansion of the particles Indeed, it shows up most clearly in the approach to 
the equilibrium, where one now expects that the plasma equilibrates the fastest m the direction 
normal to that of the field (in the color space) This is corroborated as may be seen m Fig. 
3.12 



t(fm) 

Figure 3.12: f / f eq as a function of proper time at p t = 20QMeV , £ = 0 and t c = 0.2/m for 
three different angles corresponding to cosd = 0 (solid line), cosd — .25 (dash line just below the 
solid line), and cosd = 1 (the other dash line). Note that the equlibration is fastest at 1 9 = tt/2. 


3.4 CONCLUSION 

To conclude, we have studied the production and the equilibration of a genuinely non-Abelian 
plasma with the color degree of freedom incorporated m both the source and the background 
field term m the transport equation. In the SU(2) gauge theory that we have considered here, 
the distribution function is defined m the extended phase space. We find that this approach 
recaptures in an elegant manner many of the findings of the more microscopic parton cascade 




40 


model. It has the further advantage that it indeed exhibits the color degree of freedom manifestly, 
and allows us to compute various gauge invariant quantities Significantly, we also find that the 
Abelian approximation, employed hitherto m most studies of equilibration of QGP is rather 
too drastic to be used for any quantitative analysis and comparison with the experiments The 
study also almost rules out instantaneous equilibration, and also strongly suggests that the 
collisionless limit may also not be the favoured in URHIC. 

To be sure, we have not made any comparison with the experimental findings here, for the simple 
reason that we are as yet dealing with a simpler gauge group , and ignored the gluonic component 
altogether. The indications from the present study are unmistakable, though. Indeed, the 
particle production is enhanced because of increase in the phase space available, and for the 
same reason the plasma will be cooler than the Abelian counterpart. The energy goes to the 
color degree of freedom, and does not simply heat the system as it would happen m a colorless 
plasma. If we consider the realistic 517(3) case, this feature will get further accentuated: For 
the same initial configuration, we may expect a rarer and a cooler plasma. Of course, there 
are other features which are intrinsic to 517(3): there is yet another Casimir invariant, and the 
gluon term will also have to be incorporated. They will be considered in the next chapter. 



Chapter 4 


Solution of transport equations and 
evolution of quark-gluon plasma 


4.1 Transport equations for quark-gluon plasma 

In the last chapter we have studied the evolution of a non-abelian q q plasma within the color flux- 
tube model where it was shown that the non- abelian features have a major role in the evolution 
of the system in a manner that was not captured in any of the earlier studies [18. 19, 25]. It was 
pointed out that it was unlikely that the system would equilibrate instantaneously. However, 
for simplicity, we had ignored the gluonic component and considered the simpler gauge group 
SU(2). In this chapter, we remedy both the drawbacks; we consider the gauge group SU(3) 
as is appropriate for a real plasma We also include the gluons, which have been completely 
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ignored so far We pay attention to the detailed evolution, and its prediction for different bulk 
properties of the plasma in this chapter. Signatures such as dilepton production within this 
model will be studied m the next chapter. 

There are models such as HIJING[37] and PCM[38] which study the pre-equlibrium dynamics 
via perturbative QCD (pQCD). The production and the evolution of hard and semi-hard partons 
are studied by a master rate equation in HIJING and by a transport equation in PCM. While 
they may be reliable in the study of hard partons, these perturbative approaches are admittedly 
insufficient to study the dynamics of the soft partons[37], in particular, their production. Being 
perturbative, they do not incorporate any of the non-perturbative aspects such as the formation 
of the strings and their break up, which is studied even in pp collisions. The addition of soft 
partons to hard and semi- hard partons changes the bulk properties of plasma, such as temper- 
ature and energy density. Xu et a/.[10] have observed that this addition leads to an enhanced 
suppression of which can be understood to be a consequence of increased number density 
and a lowered temperature of QGP. In our study, we employ the Schwinger mechanism[23] for 
particle production, which is quintessential^ non-perturbative. Even if one were to employ a 
perturbative version, as a time dependent electric field would require[39], which we shall propose 
in chapter six, it may be noted that an initial electric field as a classical saddle point owes its 
existence to non-perturbative processes, viz , the soft gluon exchanges that take place between 
the two nuclei. On the whole, one may expect that the pQCD based studies[37, 38] will be 
useful in the hard regime, and that the flux tube like models will be required to study the soft 
regime which will become increasingly prominent as the system expands and more and more 
secondaries are produced[40]. Indeed, it should be possible to develop a unified approach to 
study the hard and the soft components, say a la the approach of Eskola and Gyulassy[41] who 
have included the minijets in their so called chromoviscous hydrodynamics which is again based 
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on the flux tube model. This needs a separate study. 

For our purpose here, the transport equations we solve in the SU(3) gauge group to describe 
the evolution of quarks and gluons (see chapter 2 for a detail description) are: 


[p^ + Q a F*yd» + r bc Q a Aye c Q 


f(x,p,Q) = C(x,p,Q) + S(x,p.Q). 


(4.1) 


Here /(x,p, Q) is the single particle distribution function in the extended phase space of dimen- 
sion 14 which include eight color coordinates in SU(3) gauge group in addition to the usual six 
dimensional phase space of coordinates and momenta. The first term m the left hand side of 
the above equation corresponds to the usual convective flow, the second term is the non-abelian 
version of the Lorentz force term and the last term corresponds to the precession of the charge as 
described by Wong’s equation. S and C on the right hand side of the above equation correspond 
to the source and collision terms respectively. 

Note that we have to write separate equations for quarks, antiquarks and gluons since they 
belong to different representations of 517(3). 

The term gluonic source merits some explanation here. The classical background field that we 
consider here has, in contrast to the Maxwell field, self interaction. We are interested in the 
stability of the gluonic vacuum ( which is the analogue of the photon vacuum in eletrodynam- 
ics) against the fluctuations in the classical background field. An adaptation of the Schwinger 
mechanism in QED shows that the fluctuations can indeed produce the gluons, i.e. the off-shell 
classical field can spontaneously produce the on-shell radiative gluonic field (see Eqn 4.7). There 
is, therefore, no ambiguity or double counting in this process. The source term yields asymptotic 
gluonic states, and further interaction between the gluons is treated separately by a collision 
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term. 

It is very difficult to solve the transport equation written above, in general. And we have a set 
of three coupled equations here As was done in the earlier chapter we make a few assumptions. 
First of all, we admit only those potentials which can be brought to a form where the only 
surviving components are A 1111 — (A 01 , A 31 ). This choice restricts F to be “Maxwell” like, 
also pointing in the “1” direction in the color space. This restriction is not entirely arbitrary as 
it is known[42] that the non-Maxwellian configurations - where the charges, the gauge potential 
and the fields do not lie in the same direction m the group space - do not produce particle 
pairs, m general. Secondly, we require a boost invariant description of the physical quantities [2]. 
Accordingly, we demand that the distribution functions also be functions of only boost invariant 
quantities. Finally, we work within the Lorentz gauge which is implemented elegantly by the 
choice A M “ = e MI/ 3 1/ G a (r), ju, u = 0,3, with all the other components zero, r = (f 2 - z 2 ) 1/2 is 
the boost invariant proper time. 

Now we fix the magnitude of the vector charge Q a , which corresponds to the first Casimir 
invariant of SU(3). It is simply the value of the coupling constant. There is also another 
Casimir invariant, viz. d abc Q a Q b Q c , which is also conserved as the QGP evolves. There is, 
however, no way of fixing its value and the experiments presumably impose no restriction on its 
allowed values. In fact, the same holds for lattice analyses as well. For this reason, we do not 
take cognizance of this invariant, and we conveniently resolve the SU(3) charges in the polar 
coordinates: Q, — Q 0*=i sm9kCos9 t for i ^ 8, and Qs — Q TlL=i sin9k- 

As is studied by Heiselberg and Wang[43] a collision term can indeed be obtained from pQCD. 
Apart from making the equation hopelessly non-linear, this approach would be good only for 
the hard components, which are not of interest to us here. As aptly pointed by Hung and 
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Shuryak[44] recently, a microscopic description of collisions gets increasingly cumbersome and 
also unnecessary as more and more secondaries are produced. On the other hand it is admittedly 
true that there is no way to reliably obtain a collision term in the non-perturbative regime. So 
we shall employ a relaxation time approach here, where the relaxation time r c will have to be 
fixed phenomenologically, and in all probability, a posteriori. In reference[43] it is proposed 
that r c can be local and have a (weak) dependence on r. We take r c to be constant here, and 
further refinements may be incorporated after we have a better understanding of the transport 
phenomena both experimentally and theoretically. 


Within the relaxation approach, the collision term is written as 

c = zt^ALJUl. ( 4 . 2 ) 

T c 

Here f eq is the distribution function, in local equilibrium. Note that the locality can extend to 
the color space as well apart from space-time. Taking it to be an ideal gas, for simplicity, we 
write, 

2 

= exp((p^ - QM*“*K/T(r)) ± 1 (4 ' 3) 

where the +(-) sign is to be taken for fermions (bosons). Note that we allow a common equi- 
librium temperature for quarks and gluons. Here = (cosft. 77 , 0, 0, sinhrj) is the flow velocity 
and 77 is the space time rapidity given by tanhr] = z/t. With our choice of potentials, it follows 
that = 0 , so that the equilibrium distribution function can be written as: 


req 

J 9,9 


exp ((p^u M )/T(r)) ± 1 


(4.4) 


Note also that the temperature depends only on the proper time, in accordance with the Bjorken 
picture. Demanding the same of the distribution functions as well, we require that the longit- 
udinal boosts be symmetry operations on the single particle distribution. As mentioned earlier 
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the boost-invariant parameters are, apart from the color charges, 


r = (t 2 ~ r 2 ) 1 / 2 , f = (77 - y), p t = (p 2 - p 2 ) 


1/2 


(4.5) 


where y = tanh 1 (p;/po) is the momentum rapidity. 


The above set of invariant variables also serve to express the source terms for qq and gluon 
pairs which will be obtained by the Schwinger mechanism The expression for qq production 
(obtained from constant electric field) is given by[45, 46] 

27TP? 


,, , , n \ gE cos 9 1 

b q {T,t,p t ,9 1 ) = — — In 

OTT 


1 — exp 


,ct, 


) ly " 2 exp(— a£ 2 ) (4.6) 


TT 


gE cos 9i 

and for gluon pair production the corresponding term in SU(3) is given by ti j relatively stronger 
term [31, 47] 

S,(T,t,Pu» 1 ) = (3/2)S,(T,e,p,A). (4.7) 


We take g — 4 throughout our calculation. 

Now consider the third term f abc A^ a Q b -^f(x, p, Q), in the transport equation (4.1). Due to 
the restriction of the gauge potentials to the form A^ a = (A 01 , A 31 ), there is an additional 
simplification in the transport equation. In order to see that, note that equation (4.1) has the 
formal solution 


where £'(r/) is given by 


dr 1 exp(- — ) 


S'AAt'iPtiQ) [ fZKildlhSl 

Pt cosh T c 


(4.8) 


+ r dr" E(r") 

T PtT 1 Jr 1 ' 


(' = sinh -1 


(4.9) 
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Clearly, f lbc A 1 Q b -^f(x,p,Q) = 0, as S, f eq and £' only depend on Q x 

With this final simplification, we get a set of three equations, one each for quark, antiquark and 
gluon respectively. They are, explicitly, 


d 

dr 



g cos 9\E(t ) \ d 
p t cosh£ J 
f f eq 

,£ 9 |£ _ Jjhg _j_ 

T c T c 


f '< 1,3 ( ' ? £ > Pt i ^1 ) 
pt cosh £ 


(4.10) 


for quarks and gluons, and 


d 

dr 



g cos 9\E{t) 
p t cosh<f 


d 

dZ 


f q (r^,p t ,9 1 ) 


_f? , SJr J p h U } ) 
T c T c Pt cosh £ 


(4.11) 


for antiquarks. 

In the process where the field and the particles are present the conservation of energy- momentum 
is expressed in the form: 

d.TZ, + 3,77 = o, (4.12) 

where 

TZt = / PV ( 2f q + 2 /, + f g )dTd{lr, (4.13) 

and 

Tjf' = diag(£? 2 /2, E 2 / 2, E 2 / 2, -£ 2 /2). (4.14) 

Here dP = d 3 p/(27r) 3 p 0 = Ptdptd£/( 2ir) 2 , and dfl 7 is the angular integral measure in the color 
space. The factor 2 m the equation (4.13) is for two flavors of quarks. Now since energy and 
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momentum are conserved in each collision, we have: 

J p v CdFd£l 7 — 0. (4.15) 

Taking the first moment of the Boltzmann equation, integrating over the color degrees of freedom 
for f q , f q and f g and making use of equations (4.12) and (4 15), we obtain from equation (4.1) 


djf 1 ' + gE(r) J drdQ, 7 p u ^Il~JIl±M + 4 J dTdn 7 p u S q +2 J dTdtt 7 p u S g = 0 (4 16) 

In the above equation the factor 4 m third term arises because we have two separate transport 
equations for quark and antiquark, each coming with two flavors, and the factor 2 in fourth term 
is due to gg pair production, although there is only one transport equation for gluon. Putting u 
= 0 and 3 in (4 16) we get two equations, which then yield the following equation for the decay 
of the electric field: 

d -~- r dp t p] r> sinh£ r de^f, -2f q + f g } + (7r 2 /Q)d\E(r)\ 3 ' 2 = 0. (4.17) 

clt Ztc Jo Jo Jo 

Here a = a((5/2) exp(0.25/o:), a = c(g/2) 5 ^ 2 and c = Finally, C(5/2) = 1.342 is the 

Riemann zeta function. 

To solve this equation we fix the form of T(r), by demanding that the particle energy density 
differ negligibly from the equilibrium energy density, in each collision. We then relate the proper 
energy density, which is defined by 

e(r) = J dTdn 7 (p^) 2 (2f q + 2f q + f g ), (4.18) 


to the temperature by its equilibrium value, whence, 
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T{r) = i^] I/4 . (4.19) 

/ 1 

We solve the equation (4 17) numerically to study the evolution of quark-gluon plasma. 

4.2 Numerical procedure 

The numerical procedure is already discussed in the earlier chapter. We use a double self- 
consistent method. The procedure follows the scheme {T(r) ir!a /, E{r) tria {] —>{/,/, E(t)} — > 
{/, /} — y T(t) — f .. by repeated use of equations (4.8), (4.17), (4.8), (4.19), which is iterated 
until there is a convergence to the required degree of accuracy. We have in mmd the LHC 
energies, and we have taken the initial energy density as e 0 = 300 GeV/ fm 3 . Compared to 
PCM[38] where the initial particle energy density is taken to be around 1300 Gev/fm 3 at 
RHIC, our choice of initial field energy density might appear some what low However, our 
choice is guided by the estimate that the formation time for a qgp is a fraction of a fermi, as we 
explained in the previous chapter. In any case it would not be very appropriate to compare the 
two initial conditions since the ’’initiality" is only in a limited sense. Indeed, it is the energy 
density at the location of the receding color plates, and being well in the fragmentation region at 
all later times after the collision, what matters is the energy in the central region, corresponding 
to larger and larger values of r as the system evolves. It will be seen that in this region the 
results that our study yields are not unreasonable, if we make a judicious choice for the value of 
r c . We have studied, m this paper, the results for three different values of r c . For hydrodynamic 
and collisionless limits we have choosen r c = 0.001 fm and 5.0 fm. We have compared the 
results of these limitng value of r c to a realistic intermediate value r c = 0.2 fm. 
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4.3 RESULTS AND DISCUSSIONS 

The solution of the transport equations following the procedure outlined in the previous section 
allows us to determine directly the temporal evolution of the particle and energy densities, 
temperature and also the rate at which the field energy flows into the particle sector. These 
quantities are of intrinsic interest and are, in principle, amenable to experimental study via 
dilepton production. We present the results below. Also of interest are the broader questions: 
when the equilibration sets in, at what time the flow becomes hydrodynamic - with or without 
viscous flow [2, 41], etc. 

Then there are features peculiar to our model. Since the color plates are receding away from 
each other, energy is continuously pumped into the field, which subsequently decays to produce 
particles. We study the relative rates at which these two dynamical processes proceed in URHIC. 
This process cannot of course proceed indefinitely. This has already been observed by Gyulassy 
and Csernai in their study of the dynamics of fragmentation region [48]. As the plates give up 
their energy, they decelerate. The deceleration sets limits on the times up to which the model 
is valid. Indeed, the assumption of boost invariance breaks down as it is strictly exact only if 
Vpiate = C. It is safe to assume boost invariance so long as v P i ate > 0.9c[49] A simple estimate 
shows that this condition, for our choice of initial energy density, holds up to ~ 5 — 10/m for 
LHC energies. There would be other competing processes in the fragmentation region, and the 
Bjorken scenario is probably valid upto 3 — 4/m. Keeping this in mind, we have restricted 
ourselves to (proper) times < 1.5/m. Note that the results presented below will be applicable 
for RHIC energies only if r < 1/m. 

Finally, the formulation allows us to study an inherently non-abelian quantity, viz, 
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< cos 2 9 >, where 9 is the angle between the charge and the field in the group space. Note that 
its abelian counterpart = 1. Wherever possible, we have compared our results with the earlier 
517(2) study and PCM. We also display the distribution function(s) in the color space, which 
none of the other models can yield, be they flux tube based or pQCD based. All the results 
will be shown for the quarks and the gluons separately. We mention here that it is the color 
Flux-tube model within which one obtains the evolution of the mean background chromoelectric 
field, which is absent, either in HIJING or PCM. The evolution of such a background field has 
a greater impact on the acceleration of the partons present in the system, and hence on the 
observed signatures. For example the cc pair, which is produced in early collisions is acted by 
this back ground field along with cc potential, to evolve into a physical J/t/i 

4.3.1 Comparison with SU(2) results 

In the earlier chapter we had argued that the assumption of the so called abelian dominance 
which was made in a large number of transport studies[18, 19, 25] does not receive any justi- 
fication from a proper study of a non-abelian transport equation. The present study reinforces 
the same idea, as it is only to be expected. In particular, note that a determination of quantities 
such as the vacuum current or the polarization current[41] will be particularly suspect in view 
of the fact that the effective charge will only be a fraction of the true charge (see below). There 
now arises the question of the dependence of the results on N if one employs the gauge group 
SU(N) in solving equation (4.1). Although we are not in a position to make a strict comparison 
between SU(2) and S U ( 3 ) in this chapter because the SU(2) study had ignored the gluonic 
terms altogether, it is still useful to see what a limited comparison can yield. We shall restrict 
ourselves to the value r c = 0.2 fm to contrast very briefly the results of SU(2) and SU(3). 
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T(fm) 

Figure 4.2. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = .2 fm. The solid line refers to SU(3) case, and the 
dashed line to SU(2) case. 



t (fm) 


Figure 4.3: The particle number density as a function of proper time (in units of fermi), for 
r c = 2/m. The solid line refers to SU(3) case, and the dashed line to SU(2) case. 
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flow m more detail in the next section. On the whole, the plasma is denser for SU(3) due to 
increased volume m phase space. 


4.3.2 Discussion of the results 

We now discuss the results for SU(3) in some detail, for three values of r c - 0.001, 0.2 and 5 
fm. The first (last) choice corresponds to the instantaneous hydrodynamic (collisionless) case, 
m the time scale set by the initial energy density. The quark and the gluonic contributions will 
be shown separately. It is clear from Figs. 4.4 and 4.5 that the quarks and the gluons show 



x (fin) 

Figure 4.4 The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = .2 fm. The solid line refers to total energy density, 
upper dashed line to quark plus antiquark energy density, and lower dashed line to the gluon 
energy density. 

the same behaviour regarding the energy density and the number density as functions of proper 
time. The quark contribution is larger, partly because we have considered two flavors, although 
the source term[47] favors a larger rate for gluons m the color space. The situation will probably 
get reversed if we include the more correct perturbative source term for quark[39] and gluon 
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Figure 4.5. The particle number density as a function of proper time (in units of fermi), for 
r c = 2/m. The solid line refers to total number density, upper dashed line to quark plus 
antiquark number density, and lower dashed line to the gluon number density. 

production[50]. In that case, the single and three gluon production rates are of the same order 
as that of two gluons That the two sectors will continue to show the same trend may be reliably 
assumed. 

The important common feature that Figs. 4.1, 4.4 and 4.5 show is the nature of the evolution 
at times later than 1 fm. The curves suggest an approach to the hydrodynamic flow. Note 
that a hydrodynamic flow would imply that n(r), e(r) and T(t) behave like r -1 , r~ 4/3 and r~ 1/3 
respectively. We find that the corresponding exponents are —0.7, —1.23 and -0.31 respectively. 
While the temperature scaling suggests a close approach to the free flow regime, the other two 
exponents show the presence of drag[51], implying that collisions have not completely ceased. 
It may be expected that full hydrodynamic flow will take over around 2/m. 

First of all let us compare the results for quarks and gluons separately, for r c = 0.2/m before 
presenting the results for different values of r c . In Fig-4.4 we have presented the scaled particle 
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energy densites. In contrast to HIJING and PCM, we obtain lower values for energy and number 
densities for gluons than that of quarks and antiquarks together (Fig-4 5). For a complete study, 
a source term for hard parton production is also required and this term can be obtained from 
minijet production at these collider energies, following say reference^ 1]. The importance of 
the non-perturbative contribution may be gauged qualitatively by noting that the peak number 
density for quark-antiquark together(for two flavors) m this model is ~ 60 //m 3 in contrast to 
the PCM value of ~ 120 //m 3 (for three massless flavors). The energy density for quark sectors 
in PCM is also larger by the same factor, implying that the average energy per particle is of the 
same order in both the cases. 

Let us now consider < cos 2 9 > (where 9 is the angle between color charge and the chromoelectric 
field m color space). This quantity is a good bench mark to characterise the 'non-abelan ness of 
the system, and we present our results in Fig-4 6. This quantity is gauge invariant and physical. 



i(fm) 

Figure 4 6: < cos 2 9 > as a function of proper time (in fermi) at r c = .2/m for quark (solid line) 
and gluon (dashed line) 

It may be seen that this value saturates at 9 ~ rr/4, for both quarks and gluons. This value was 
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always unity m an abelian plasma, as 9 = 0 in that case. 

This has a direct impact on the equilibration of the plasma. The equilibration is faster around 
9 ~ tt/ 2 where the background field effect is zero, and is slower at 9 ~ 0 This is under- 
stood as follows. At larger value of 9, say around tt/2 , there is no acceleration of partons by 
the background electric field. Only collisions are present at this angle, and hence the rate of 
equilibration is faster. On the otherhand when the angle # = 0, the acceleration of partons by 
background field retards the equilibration of the partons. At any intermediate values of 9 , the 
rate of equilibration lies m between these two extreme values. So this average value, which is 
purely due to the non-abelian nature of the theory, has a major role in the equilibration of the 
plasma. This is more clearly shown m the distribution functions which carry the color degrees 
of freedom explicitly(Fig-4.7 and Fig-4.8). 



Figure 4 7; f a /f‘ q as a function of proper time at p t = 300 MeV, £ = 0 and r c = 0.2/m, for 
different values of 9. Solid line corresponds to 6 = tt/2. 


We conclude the discussion at r c = 0.2/m by making an interesting observation. As can be seen 


m Fig-4.6, the value of < c.os 2 9 > fluctutates around its mean value for both quarks and gluons. 
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x (fm) 

Figure 4.8: f 3 jf tq as a function of proper time at p t = 300MeV, £ = 0 and r c = 0.2/m. for 
different values of 6 Solid line corresponds to 9 — tt/2. 

with similar fluctuations in other quantities as well - although it is not that pronounced m them. 
While it certainly indicates that the distribution is approaching its equilibrium value asymp- 
totically, the question is whether the fluctuations would persist even as the system hadronises. 
Unfortunately, it is difficult to make any definite assertion at this stage. Indeed, as pointed 
earlier, evolution of the system beyond 3 — 4 fm requires a full treatment that goes beyond the 
assumption of the Bjorken scenario in the central region. It is also not known when the hadron- 
ization exactly sets in. However, if one assumes that the fluctuations do indeed persist, it may 
quite well happen that it will manifest as an anisotropy in the parton momentum distribution. 
This has been pointed out in a study of the related ‘color filamentation’ by Mrowczynski[52] 
It is also possible[53], although we do not know in what manner, that these fluctuations in 
the hadronic phase show up as disoriented chiral condensates. It is only a more complete and 
rigorous study that can settle the status of these speculations. 


Now we present the data for the other two values of r c , namely r c = 0.001/m and t c = 5/m, 
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corresponding to hydrodynamic and collisionless limits respectively. These results are then 
compared with ~ c = 0.2/m. Consider first the number density n(r). This quantity is seen to be 
sensitive to the values of r c . As can be seen from Fig-4 9, n(r) is larger m the hydrodynamic 



Figure 4.9: The particle number density as a function of proper time (in units of fermi), for 
t c = 2/m (solid line), for r c = 0 001 fm (upper dashed line) and for r c = 5/m (lower dashed 
line). 

limit, where the maximum value is around 120 per fm 3 , at r ~ 1.0/m. As r c is increased to 
0.2 fm we find a lesser value, the maximum value being around 80//m 3 . This trend continues 
for collosionless limit where particle density is still less. The increase in number density as r c 
is decreased is expected in the context of classical theory. This is because, the collision time 
m any classical non-equilibrium theory depends roughly on the inverse of the number density, 
apart from the other factors. The behaviour of number densities in Fig-4.9 for different values 
of r c in oui; calculation also reflects the above fact. 

In Fig-4.10 we present the scaled energy densities for different values of r c . Unlike the number 
density, there is no strong dependence on r c . This means that the average energy per particle is 
different for different values of r e . We find that for r c = 0.001/m, the maximum average energy 
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t (fm) 

Figure 4.10: The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for r c = .2 fm (solid line), for r c = 0 001/m (upper dashed 
line) and for r c — 5/m (lower dashed line). 

per particle is around 1.0 GeV, whereas it is around 2.0 GeV for r c = 0.2/m and 4.0 GeV for 
r c = 5.0/m. Note that these high energy deconfined partons can produce secondary partons, 
such as strange quarks and also can break a fully formed J /ip as analysed by short-distance 
QCD in reference[54]. 

Let us consider the electric field. It may be seen from Fig-4.11 that the decay of the field is 
very slow at r c = 0.001/m, compared to r c = 0.2 and 5 fm. Only 30 percent of the field 
has decayed m the hydrodynamic limit in comparison to collisionless and intermediate limits, 
where the decay is very large. At the other end in the collisionsless limit, the decay is at a 
slightly slower rate in comparison to r c = 0.2 fm, implying that the maximal conversion rate 
is for t c around 0.2 fm. As mentioned earlier, one peculiarity of the model is that the field 
energv is continuously created due to the recession of the plates even as the field itself decays 
to produce particles. To study this, in Fig-4.12 we present field energy per unit transverse area 
as a function of ordinary time. For r c = 0.2 fm this energy is much less than that at t c = 5 
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Figure 4.11: Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for t c = .2/m (solid line), for r c = . 001/m (upper dashed line) and for r c = 5/m (lower dashed 
line). 


fm and 0.001 fm. This demonstrates that the conversion is indeed dominant at r c = 0.2 fm. 
This is more clearly displayed in Fig-4 13, where we have plotted the ratio of particle energy 
per unit transverse area to field energy per unit transverse area. 


In Fig-4.14 we present the evolution of temperature at r c = 0.2 fm. The maximum temperature 
we obtain in the color flux-tube model is around 300 MeV. Finally, a brief comment on the 
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Figure 4.12: The field energy/unit transverse area(m units of GeV//m 2 ) for r c = 0.2 fm (solid 
line), for r c = 0.001 fm (upper dashed line) and for r c = 5 fm (lower dashed line), as a function 
of ordinary time (in fermi). 



Figure 4.13: The ratio of particle energy per unit transverse area to field energy per unit 
transverse area for r c = 0.2 fm (solid line), for t c = 5 fm (upper dashed line) and for r c = .001 
fm (lower dashed line), as a function of ordinary time (in fermi). 
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distribution functions for quarks and gluons with explicit color dependence (Fig-4.7 and 4 8). 
As discussed earlier, it may be seen that the equilibration is fastest for 6 = tt/2. where there is 
no background effect. At 6 = tt/2 the thermal equilibration occurs at a common time r ~ 1 0 
fm for both quarks and gluons. 


4.4 Conclusion 

We have studied the production and the equilibration of a genuinely non-Abelian plasma with 
the color degree of freedom incorporated in both the source term and the background term 
m the transport equation. With the realistic gauge group SU(3) that we have considerd here, 
the distribution functions for quarks and gluons get defined in the extended phase space of 
dimension 14 We have added the gluonic component into the color flux-tube model, which was 
so far absent 

The role of the color degree and that of gluons was found to be substantial, and non-trivial. The 
plasma is denser than the SU(2) counter part but is comparable to the Abelian case. However, 
it is much cooler than the abelian plasma, where the corresponding temperature is ~ 800 MeV. 
The value of the effective charge is larger than that of its counterpart in SU(2), contrary to naive 
expectations. The peak energy per particle is ~ 2 GeV which indeed is the demarcating scale[37] 
between soft and hard processes. In short, the non-perturbative aspects of the evolution of QGP 
offer a rich variety of results which will have to be combined with the perturbative studies in 
order to abtain a complete description of the production and evolution of quark-gluon plasma. 



Chapter 5 


Dilepton production from a thermally 
equilibrating quark-gluon plasma 


As mentioned earler, some of the major proposed signatures to detect quark-gluon plasma 
in URHIC are 1) Jjt' supression[5], 2) electromagnetic probes such as dilepton and direct 
photon production(6, 7], and 3) strangeness enhancement [8] . It has been suggested that for an 
equilibrated quark-gluon plasma, J/xp supression, due to screening, is a good probe[5]. However, 
due to the uncertainties in the study of production and evolution of QGP m these is experiments 
the results obtained by screening can not be taken too seriously. There are also calculations 
of J/tJj suppression in the equilibrating quark-gluon plasma[9, 10] using short-distance QCD 
which are different for different models (see chapter 7). The supression of J/ip. infact, has been 
observed in reactions where there is no quark-gluon plasma phase, such as in p-A collisions and 
in collisions of light nucleus[55]. This suppression is well explained by the nuclear absorption 
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of J/xp[ 56] Although, recently, NA50 collaboration[57] reports an excess in the supression of 
J/0, it is still not clear if a quark-gluon plasma phase has formed in this collisions. There are 
proposals that the data are explained by a deconfined partonic medium[58], or by a medium with 
high density[59]; but there are also other calculations which explain the data without assuming 
any quark-gluon plasma phase[60]. As far as the J/xb supression in nucleus-nucleus collision is 
concerned, many aspects of it has to be studied in greater detail before unambigious conclusions 
can be drawn about the existence of quark-gluon plasma. It thus appears that there is a need 
to study more than one signature if one has to detect QGP. With that in mind we investigate 
another signature, i.e. dilepton emission. 

Dileptons and single photons have long been proposed as useful probes of the plasma[7]. Once 
produced, they hardly interact with the strongly interacting matter and thus carry the details of 
the circumstances of their production. There has been, however, uncertainties m the calculation 
in their spectra. This is due to lack of detailed knowledge about the space time evolution 
of quark-gluon plasma. Recall that the various stages by which the complete evolution of 
quark-gluon plasma is described in URHIC are, i) pre-equilibrium, ii) equilibrium, where one 
actually studies the equilibrated quark-gluon plasma, iii) cooling and iv) hadronisation. The 
pre-equilibrium stage of the collision which leads to thermal and then chemical equilibrium has 
a crucial role to play. From this point of view it is necessary to study what happens to the 
dilepton emissions from different stages of evolution of QGP, rather than estimating it in an 
equilibrated quark-gluon plasma. 

We study the dilepton production in an equilibrating QGP employing the distribution functions 
determined in the previous chapter. The dominant process which produces a dilepton pair l + l~ 
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is 


q + q 7’ ->• / + + / 


(5.1) 


where 7* is the intermediate virtual photon. The dilepton emission rate in a space time volume 
d A x is written as 


dN 

d 4 xd 4 P 


1 

(2? r) 6 


J d 3 p 1 d 3 p 2 dQf(x,p u Q)f(x,p 2 ,Q)v re icr(M 2 )S A (P - pi - p 2 ). 


(5 2) 


Here is the four momentum of the lepton pair, with P? its transverse momentum and Mr(= 
\j M 2 + Pt 2 ) its transverse mass. M is the invariant mass of the dilepton pair (M 2 = P^P^) 
and v re i = M 2 /2 EiE 2 is the relative flux velocity of quark and antiquark pair in the above 
process. The distribution function of quark and antiquark is a function of the color charge Q. 
dQ is the integral of the the color part of the phase space. The dilepton production cross section 
for the above reaction is 
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where F q = N 2 jjJE s t) and m t is the mass of the le P ton - N 1S tiie color averaging factor 
(corresponds to the volume of the compact (color) space), N s is the spin degeneracy ( N s = 2s + l) 
and ej is the fractional charge of the flavour. For a dilepton pair of large mass the emission rate 
in the midrapidity region (Y = 0) of lepton pair is given by [61]: 
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withp t2 = - UhMchtn ~ 6) + P?i» sh & = j^{M T shr]-pnsh^) and p ± = \M 2 [M T ch{p- 

T p T \~ 1 . In equation (5.4) we have used d A x = ■nPPdrTdq (where R is the radius of the 
radiating region). a(= 1/137) is the coupling constant of the electromagnetic interaction. 
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It has been shown[62] that the thermal dilepton rate does not depend on M and Pt separately, 
but only on Mj. This is known as Mt scaling in QGP. However this Mt scaling is violated 
m an equilibiating plasma, as can be seen from equation (5.4). In this case the dilepton rate 
depends on both Mt and Pt This would also be the case if one takes the transverse expansion 
of the plasma into account [63]. 

In our calculation we have taken R = 7/m corresponding to Au-Au collisions. We have con- 
sidered two different cases corresponding to two different relaxation times, r c = 5/m and 
r c = 0 2/m. The relaxation time r c = 5/m corresponds to collisionless limit and r c = 0.2/m 
corresponds to a more realistic limit of equilibration[38]. We have calculated the dilepton rate 
as a function of Mt for Pt = 0 and 1 GeV respectively. In Fig-5.1 we have presented the 
dilepton rate at r c = 0.2/m. As’ can be seen from Fig. 5.1, the dilepton yield becomes smaller 



Vt^GeV) 

Figure 5.1: Dilepton rate from pre-equilibrium stage, as a function of Mt for r c =0.2 fm. 

for higher values of transverse momenta. In Fig-5.2 we have shown our results for r c = 5/m 
which corresponds to the collisionless limit. The rate at higher Mt is found to be larger than 
that at t c = 0.2/m. This is because the average energy per parton is higher at r c = 5/m than at 
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t c — 0.2/m as observed in the last chapter. This is because in the collisionless limit the energy 
per parton is larger (~ 4 GeV) than that at r c = 0.2 fm (~ 2 GeV). 

In Fig-5 3 and Fig-5.4 we have compared our results with the corresponding Drell-Yan predic- 
tions. The Drell-Yan results are taken from the reference[28]. It can be seen from Fig. 5.3 that 
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Figure 5.2: Dilepton rate from pre-equilibrium stage, as a function of Mt for t c =5.0 fm 
(collisionless limit). 

the dilepton rate from pre-equilibrium stage is larger than the Drell-Yan production for very 
small transverse momentum (Pj — 0 GeV). This is true when Mj < 2GeV. For Mj > 2 GeV, 
the Drell-Yan production dominates over the pre-equilibrium dilepton production. For a large 
transverse momentum ( Pf — \G& V") of dilepton pair, the production from the pre-equilibrium 
phase is lower than the Drell-Yan production in the whole range of dilepton transverse mass. 
This can be seen from Fig-5.4. This is contrary to the earlier findings [64-66] where it is shown 
that dilepton production from the pre-equilibrium stage dominates over the Drell-Yan produc- 
tion. However the production of partons are not too hard in this model, and one expects to 
have such a low rate in the dilepton spectra than the Drell-Yan emission. 
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Figure 5.3. Drell-Yan and pre-equilibrium dilepton rate as a function of Mt for r c =0 2 fm 
(P T = 0 GeV). 


Summarizing the chapter, we have calculated the dilepton spectra in ultra relativistic heavy-ion 


collisions within color flux-tube model, with non-abelian features explicitly incorporated. In 


contrast to the earlier studies, we have found a smaller rate in dilepton production in the pre- 


equilibrium stage at high Mt and Pp than the Drell-Yan production. The production is larger 


in the collisionless limit of the plasma. An experimental measurement of the dilepton spectra m 


URHIC can shed light on the collision time (r c ) and possibly on the initial field energy density 


that we have assumed in our model. 






Chapter 6 


A new source term for quark production 
from a space-time dependent 
chromoelectric field 


6.1 Introduction 

Over years there has been several devlopement on the production of charged particles from 
electric field, a phenomenon which was discovered nearly forty years ago by Schwinger [23]. By 
now the production of particles from a constant background electric field is extensively studied 
both theoretically and experimentally[67]. However, the subject matter of particle creation 
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from a time dependent field, and the dynamics of the ffield-particle' system is relatively new 
and is not fully solved. This aspect has been a subject of renewed and central interest in the 
context of heavv-ion collisions and QCD due to the suggestion that the receding nuceli might 
pioduce a strong chromoelectric field which creates quark/ anti-quark pairs and gluons[19]. As 
the particles (quark/anti-quarks and gluons) are continously produced from the background 
field, the field becomes time dependent during the evolution of the system after the two nuclei 
collide with each other. Hence the source term for the production of particles which is taken in 
chapter .3 and 4 incorporating a constant field in its calculation is not strictly applicable and a 
determination of the particle production from a time dependent electric field becomes necessary. 
Indeed a complete quantum mechanical treatment of the problem is rather too involved. For our 
purpose, it should be sufficient to derive how a space-time dependent field produce particles as 
a function of time. The dynamics of the field is of course is coupled to that of the particles, and 
may be determined in a consistent (although not complete) manner by solving the transport 
equations as we have done in the previous chapters. 

In this chapter we derive a source term for quark/anti-quark and gluon production from a time 
dependent electric field which is needed in the non-abelian transport equation[45], with the 
purpose of eventually studying the production and equilibration of quark-gluon plasma in ultra 
relativistic heavy-ion collisions. As described in the preceeding chapters the source term one 
needs is, dt J^p t , which is the probability for particle production per unit time per unit volume 
per unit transverse momentum in the phase space. 

Before deriving a source term for a time dependent electric field there are several physical 
features that needs to be clarified. First of all Schwingers’s mechanism[23] incorporating a 
constant field in its calculation is an exact one loop non-paeturbative result, which can be 
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understood m terms of a semiclassical tunneling process[6S]. However, for a time dependent 
electric field such a calculation is not necessary and particle production can directly take place 
by the perturbative mechanism. It can be further noted that a constant field can not produce 
pairs by this perturbative mechanism. Secondly, it is not possible to obtain a source term 
(dic^xePpf) m usua l phase space from a quantum processes This is because any attempt to 
obtain such a source term within zero uncertainty will violate quantum mechanics. Although a 
source term m the usual phase space has been proposed[68] by rewriting the Schwinger’s source 
term as an integral over pf , such proposals are ad hoc and are inconsistent with the minimum 
uncertainty requirement. In this chapter we overcome this problem by working in the coherent 
basis of coordinate, momenutm and color. 

A coherent state (|a >) is the eigen state of anihilation operator a(x,p) defined by [69] 

a\a >= a|a > (6.1) 

with a being its eigen value the real part of which is proportional to the expectation value 
of cordinate operator and the imaginary part is proportional to the expectation value of the 
momentum operator. W.r.t this coherent state the expectation value of coordinate operator x 
and momentum operator p can be written as 

X =< a\x\a > (6-2) 

P =< a\p\a > (6.3) 

with AXAP = h j 4?r , which satishfy the minimum uncertainity. Here (AX) 2 = < a\x\a > 2 — < 
o|i 2 |a > and (AH) 2 = < a|p|a > 2 — < a|p 2 |a >. We have used this coherent basis Jq > (= 
\X, P >) to derive a new source term which does not violate principles of quantum mechanics. 
To derive a source term for parton production from a time dependent electric field and to further 
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study its dynamics we have restricted our calculation to a finite time interval in the S matrix 
expansion (see below). 

In the next section we derive the source term for the case of a purely time dependent electric 
field. Pending the determination of E(t) from a self consistent study, we have employed a specific 
form for the electric field to get an understanding of the dynamics of the particle production. We 
have compared the results of this source term with that of Schwinger’s source term (obtained 
for a constant field) by putmg the same time dependence in the field (the time dependence m the 
electric field in Schwinger’s source term was put by hand m chapter 3 and 4). In section 3 of this 
chapter we have derived a source term for quark production from a r(= yW— z 2 ) dependent 
elctric field. These source terms will be used in the transport equation to study production and 
equilibration of quark-gluon plasma in ultrarelativistic heavy-ion collisions. 


6.2 Source term for a purely time dependent electric field 

We start our derivation with a determination of the Feynmann amplitude for the production of 
quark from a time dependent electric chromoelectric field. The Feynmann amplitude is written 
either in momentum space or coordinate space. However, we need a source term ( ) i n t ^ ie 
usual phase space of coordinate and momentum from this amplitude. For this purpose, we have 
projected the Feynmann amplitude (m momentum space) to a coherent basis of coordinate and 
momentum. As described earlier this coherent state is the eigen state of anihilation operator 
and the expectation value of the coordinate and momentum with this coherent states satishfy 


the minimum uncertainty. 
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For a particle at mean position X, having mean momentum P, and color charge Q the amplitude 
in the coherent basis can be written as 

<Q,X,P\S\Q>= I d z k <Q,X,P\kxk\S\0>, (6.4) 

where |A', P,Q > is the coherent state of cordinate, momentum and color. Here the coherent 
state |Q > for the color is defined as | Q >= exp(iQ a T a ) |0 > (see Eqn 6.14). To the leading 
order in perturbation series we have 

$ ~ J dt'd z xLx(t',x). (6.5) 

with C%(t , x) being the interaction lagrangian density. The time (t = 0) refers to the initial 
state where there is no production (this replaces t = — oo in the usual scattering processes). 
For our dynamics here we have limitted the time integral to a finite limit t and not to co. For 
quark/anti-quark pair production the interaction lagrangian density is 

Lx = (6 6) 

with i /;, being the Dirac field for quark and A a being the external color field. In the above 
equation T* are the generators with a = l..,.,8 and z, j = 1,2,3. For quark/anti-quark pair 
production the corresponding amplitude is 

< Q\i Oh, -Yi , Ah, Pi- P>\S\Q >= J d z kid z ki < Qi, Q 2 , Xi, JY 2 , Pi, i^l^i, J > 

< z,j,Si.s 2 ,ki,k 2 \S\0 >, (6.7) 


where ,Y t , AQ are coordinates P x , P 2 momenta and Si and s 2 are the spin indices of the quark 
and anti-quark respectively. Here < i,j,Si,s 2 ,ki, ^I^IO > is the usual Feynmann amplitude for 
quark/anti-quark pair production with i and j being the color index of quark and anti-quark. 
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After summing over spin and color the absolute square of this amplitude is written as: 


|< QuQ2,X u X 2 ,P u P 2 \S\Q >| 2 = J d%d 3 k[d%d% < Q u Q 2 ,X u X 2 ,P u P 3 \k u k 2 l i,'', 

< l,m,k l ,k 2 \Q l ,Q 2 i Xi,X 2 ,Pi,P2 >*< 0|5|fci, k' 2 , s u s 2 , /, m >*< z, y, si, s 2 , ki^l-FlO > 


In our calculation, we ha\e made a suitable gauge choice such that = 0. With this 

choice the potential .4 M in the the above equation is related to the electric field by the simple 
expression. 


E*(t) = - 


dA^jt) 

dt 


(6.9) 


Substituting 

A(t, x) = J d 3 KA(t , K) exp {—iK ■ x) 


( 6 . 10 ) 


we obtain (for massless quark) 


— 'u < 0|*5 |^i i k 2 > < ky , > — f(t\, ki, k 2 , A: 2 ), (6.11) 

where f{t[,t 2 ,ki,k 2 , k [ , k 2 ) is given in the appendix for a fixed value of a in A M “. The expression 
foi color sum (C(Q)) is not included in f(ti,t 2 ,ki,k 2 ,k[,k'2) and is presented separately (see 
Eqn (6.14)). With the above fourier transformation (Eqn (6.10)) the time dependence in A 
appears finally in the source term. In most of the earlier studies[71] the source term is derived 
with some special form of the electric field. Such source terms are not necessarily applicable 
m the transport equations to study the production and equilibtation of quark-gluon plasma. In 
our study, here no specific form for electric field is assumed. 

For the case of a single quark we project Fevnmann amplitude to a coherent basis \X, P,Q > 
rather than projecting it to | X u X 2 , P u P 2 , Q l, Q 2 >■ The later is needed if we are interested in 
pair (quark/anti-quark) production and corelations. For quark production at mean position X 
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having mean momentum P and color charge Q we have 


|< Q, X, P|5|0 >| 2 = 


k <Q,X,P,\k u : ixj k[\Q,X,P>* 


< 0\S\k[,k,j, l >* < l,i,k u k\S\0 > . ( 6 . 12 ) 


The above equation can be written as 


|< Q,.Y, P|A |0 >|' 


k<Q\i><j\Q >T t a iT' b h < X , P, \k\ > < k[\X,P > 


< /n,£|P a |0 > < 0|5 6 |fci,yfc >(*6.13) 


As was done m the previous chapters we fix the direction of the electric field and we take it to 


be along 8 direction in the color space (a =8 and b= 8 ). The color factor C(Q) becomes, 

C(Q) = S tJl , < <51* > < J\Q >*T$T* % = ^[1 + 3| < Q|3 > | 2 ], (6.14) 


with |Q > = ex/i(iQ a T a )|0 > 


Following the reference[70) we arrive at the expression < Q |3 > = cos 6 ex 5 ( 17 ); so that C(Q) = 
^[1 + 3 cas J /?]. Tins is the univerasl form of the color dependence independent of the nature of 
I hr electric field. Here 0 < 0 < tt/2 and 0 < 7 < tt. 

For a purely time dependent potential A(t.ki + fc 2 = 0) we obtain 

|< Q,A',P|S|0 >| 2 = C(Q) I d 3 M 3 fcl2d 3 fc < X,P, \h > < k[\X,P> ’ 

[ dt 1 f dt 2 f{t l ,t 2 ,k l ,k,k[,k)exp(i(k° J rk 0 )t l )exp(-i(k[ 0 + k 0 )t 2 )5 3 (ki + k)5 3 (k[ + kf.Q.15) 
J 0 Jo 

The coherent state wave function < X , P\k > which satisfies minimum uncertainty is 


< X, P\k >= N exp(—(P — k) 2 /2Q 2 ) exp(ik ■ X), 


(6.16) 
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with N being the normalisation constant. Putting the the above form for the coherent state m 
equation (6.15) we get 


!< VPIS10 >| ! = N*C(Q) J cPkexp(—(k — P) 2 /Q 2 ) f it, £ dt 2 

expC2 1 fc°(i 1 - i 2 ))-J_^[P(A( (l ) . A(i 2 )) - (fc A(t 2 ))(t-A(t 2 ))], (6 17) 

with k — | A. | (for massless quark) Without any loss of generality one can assume z-axis to be 
along the direction of A and obtains 


|< Q, A, P|/>|0 >| — C[Q)N 2 J d 3 kexp( — (k — P) 2 /Q 2 ) dti J dt 2 (A(ti)A(t 2 ))exp( 2 i\k\(ti — ti) 


-[(1 - cos 2 9 ^ | 


is: 


wit.h <9* being the angle between z-axis and k. This is the probability of quark production per 


unit volume per unit momentum having a color charge Q from a time dependent electric field 

E[t) related to A(t) by Eqn ( 6 ^J t The required source term is then 

,s 


dN 


dtd 3 


(iV _ r 

XdPPr ~ J 


dP, 


d\ < Q,X,P\S\0 > 


dt 


(6.19) 


To make sur<' that a constant potential and a constant field (which is obtained from a linear time 
dependent, potential (Eqn. (6.9)) do not contribute to the above time dependent source term 
we project out only the time dependent part by the replacement A(t) A(t) — A(0) — tA'(O) 
in the Kipi. (6.19). After singling out the contributions from this constant potential and the 
constant held we obtain a genuinely time dependent source term from Eqn. (6.19). 

Now by putting some specific form of the electric field we compare the results obtained by this 
source term with that of Schwinger’s source term. For the choice A(t) — exp(— t/to) we present 
the results of both the cases below. We have choosen Qx a,ud Q 3 m the coherent state wave 
function (Eqn. 6.16) corresponding to a 1T1 dimensional expanding plasma. For Qx the scale 



80 


is the transverse diameter of the nucleus (~ 14 fm for Au nucleus) and for Q 3 it is determined 
from the total length trailed by the two nuclei (2z or 2t, as nuclei travel almost at a speed of 
light,) In Fig-6.1 we have plotted the production rate for P T 0-5 GeV as a function 

of t The production rate increases as time increases in contrast to the Schwinger’s results 



Figure 6.1: The rate of quark production {GeV 2 )) obtained by the new source term 

from a time dependent electric field for Pt = 0.5 GeV, f 0 = 0.5 fm and 9 = 0 

which is shown in Fig-6.2 for the same input parameters. This is because, in the newsource 
term we have taken a finite limit in the time integration in the feynmann amplitude. As time 
increases, the absolute value of this amplitude increases and hence a large production rate is 
obtained. However, the schwinger formula is valid for a constant field. But for a comparison 
with our new source term we have taken the same time dependence in the electric field in both 
the case's. In Fig-6.3 we present the results for Pt = 1.0 GeV. The production rate decreases 
as P r increases, which is in accordance with the Schwinger’s result (presented in Fig-6.4). 
This is more clearely seen m Fig-6.5 where where we have plotted the production rate of the 
new source term for Pt = 5.0 GeV. To get a feeling on the total production of particles rather 
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Finnic* 6 2: The rate of quark production ( dtd3 (<j e V 2 )) obtained from Schwinger’s source 
trim with the same time dependent electric field (as described in the text) for Pj = 0.5 GeV, 
/o - 0.5 f'm and 0 = 0 

than its production rate we have plotted the total production ( d i£#p- ) Fig-6.6 and 6.7 for 

p r - u.5 and 5.0 (le.V up to t = 2.5 fm. It is seen that total number of particle production per 
unit volume per unit transverse momentum is an increasing function of time. This has a strong 
impact, on the particle production in relativistic heavy-ion collisions where the plasma would be 
very ilen.se even at a later time (up to 2.0 fm). Such a calculation of the quark-gluon plasma 
production and its evolution within the frame work of non-Abelian transport equation will be 
performed and reported elsewhere. 
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P igur(“ 0..5: I lie rate of quark production ( "A uFTiP p r ( Ge V 2 ) ) obtained by the new source term 
from a tune dependent electric field for P T = 1.0 GeV, t 0 = 0.5 fm and 9 = 0 



Figure 6.4: The rate of quark production (^f ^(GeV 2 )) obtained from Schwinger’s source 
term with the same time dependent electric field (as described in the text) for P? = 1.0 GeV, 
h j = 0.5 fm and 0 = 0 
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Figuie fi 5 Fix* rat<> of quark production ij^i^rp^iGeV 2 )) obtained by the new source term 
fiorn a time dependent, electric field for P T = 5 0 GeV, t 0 = 0.5 fm and 9 — 0 



Figure tUi: I'he quark production ( d 3 T&p^ (GeV )) obtained by the new source term from a 
time dependent electric field for Pt == 0-5 GeV, to = 0-5 fni and 9 — 0 
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Figure 6.7: The quark production p- (GeV )) obtained by the new source term from a 

time dependent electric field for Pr = 5-0 GeV, to = 0.5 fm and 9 — 0 
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6.3 Source term for a r dependent electric field 


As seen m the previous chapters, for a 1+1 dimensional expanding quark-gluon plasma we 
have employed Bjorken’s boost invariance formalism according to which the longitudinal boost 
are the symmetry operations on the particle production. The boost invariant parameters are 
r = sj{t 2 — z 2 ) y Pj and f = 77 — y with 7 7 = tanh~ 1 (z/t) and y = tanh~ l {P 2 ! Po). The source 
term for particle production ( d $x% p^ -(r, £, Pj)), needed m the transport equation depends only 
on these boost invariant parameters. In this section we derive such a source term which is 
needed in the transport equation to study production and equilibration of quark-gluon plasma 
m ultrarelativistic heavv-ion collisions. 


In terms of these variables we get dtdz = r dr drj and 

S = f t' dr' dr)dxdyCx{r' , 17, x, y). 

Jo 


(6 20 ) 


( r 


(this equation is analogous to Eqn. {?) for a purely time dependent electric field). To be 
consistent with the previous chapters we work in the Lorentz gauge A 1 = 0 with A 1 = 




d„G(r) 


where /*, u = 0 , 3 . Due to this choice the chromo-electric field depends only on the 


boost invariant quantity r which is given by 


e{t) = {AA-Ug(t) = [J-A]g\t). 


Following the same procedure as was 


dr 2 ' rdr J ~ w 1 dr r 

done for the time dependent case above, we get 


( 6 . 21 ) 


\<Q,X,P\S\Q> 


| 2 = J d%d%d z k <Q\%X3\Q >TZr h lj < X , P, \h > < k[\X, P > 

[ T dr x n f T dTrtdmdmgih, k, k[, k, n, r 2 , m , V 2) exp(*(*? + k^ncosh^) 

Jo Jo 

exp (-i(fci° + k 0 )r 2 cos/u] 2 )exp(-i(A: lz + k^sinh^exviKKz + K)T 2 sinhr, 2 ) 

S 2 ( Jfefr + k T )5 2 (k[ T + %22) 
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whete the expression for g(k x ,k, k[, k, r 1; r 2 , 771, rj 2 ) is given in the appendix for massless quark. In 
the above equation we have used the fourier transformation A{t , z, X t ) = f cPK t A(t, z, K t ) exp(-i K t - 
Xt), with A(t, z, Aj = 0) for the vector potential which depends only on t and z. Eqn (6.22) 
can be further simplified and is written as 

|< Q,X,P\S\0 >| 2 = N 2 C{Q ) J dk lz dk[ z d 2 k T dk z ex V {-{k T - P T ) 2 /Q 2 T ) 
exp(-((k lz - P 3 ) 2 + [k u - P3) 2 ) / 2 Ql) dr X T x J dr 2 r 2 dr]idrj 2 

g{k u k, k[, k,T U T 2 ,rh,ri 2 ) exp(i(k° + k°)TiCoshr)i) exp(-i(k[° + k°)T 2 coshrj 2 ) 

exp(— i(k Xz + k z )T i sxnhr] 1 )exp(i(k[ 2 + k z )r 2 sinhi 72) exp(— i(k u — k[ z )X 3 ), (6.23) 

with k° = \J(k? + yfc 2 ), k\ = \/{k^+~kf z ) and k'° = + kfi.). Substituting X3 — rsmhr/ and 

P3 = Pjsinhy with £ = 77 — y we get 

|< Q,X,P\S\0 >| 2 = N 2 C(Q) J dk lz dk' u d 2 k T dk z ex V {-{k T - P t ) 2 /Qt) 

exp( — ((*1, - Prsinhy) 2 + (k lz - P T smhy) 2 )/ 2 Ql) f dr x r x f dr 2 r 2 dr] X dr] 2 

Jo J 0 

S{k u k, k[,k, g u g 2 )G l (T l )G'(r 2 ) exp(i(&? + fc°)r 1 cos/i77 1 )exp(-7(ki° + k Q )r 2 coshr ] 2 ) 

exp {-Jk lz + k z )Tismhr}i)ex'p(i(k' lz + jfe,)r 2 a*n/»772)exp(-i(A:u ~ k' lz )Tsinh{ £ + y)) (6 24) 

where 

S(k x , k, k[, k, 771 , 772 )G^ (ti )6? / (t 2 ) = g(ki, k , Aq, k, Ti, t 2 , 771, 772). (6.25) 

Integrating this expression over P 3 and taking derivative w.r.t. r we arrive at the new source 
trem, 

inriTv ' = Pt t f d v cosh v l< l 5 l° >l 2 - ( 6 - 26 ) 

dtd 3 xd 2 Pr dr J 

which depends only on the boost invariant quantities r, £ and Pr- This source term will be 
used in the transport equation to study the production and equilibration of quark-gluon plasma 
m ultrarelativistic heavy-ion collisions. This work will be taken up separately. 
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6.4 APPENDIX 


For m. = 0 one obtains 


f(ti, t 2 , ki, k 2 . k[, k 2 ) — — (/i + /2) 

JD 


(6 27) 


where 


B = ^/k°k^k°k[° 


(6.28) 


and 


/i = pi • A)(& 2 • A‘)(4k[°k' 2 ° - (k[ k 2 )) + (k[ • A)(fc 2 • A*)(A: 1 • fc') + (h A)(k[ • A m )(k[ • fc 2 ) 
-(k[-A)(k' 2 A“)(ki ■ k 2 )] + pi • A*){k 2 A)(4k' 1 °k 2 0 — (k[ ■ k 2 )) + (k[ • A*)(k 2 • A)(k\ k 2 ) + 
(ki ■ A*)(k[ ■ A)(k[ • k 2 ) — (fci • A*)(k 2 ■ A)(ky ■ k 2 )] + (A • A*)[4A;( k 2 (k± ■ k 2 ) + ( ki ■ k x ){k 2 ■ k 2 ) 


+(*, ■ *,)(*; ■ K) - (fc, ■ k' 2 )(k[ ■ k 2 )j + [((*! • aka; ■ a-) - (t, • A-)(t; - ajx*, • k' 2 ) + «k 2 


(*4 ■ A) - (fcj ■ A)(4 . A*))(fc x *405-29) 


and 


I 2 = (fc, • A)(Jfc 2 • A-)[2i-;°t' 0 - (k[ ■ k ' 2 ) ] + (*,, A)(*4 ■ A-)[(*2 • *4) - 2*4 V] + (*9 • A”)(ki • A) 
[2k[%° - (k[ - k' 2 )\ + (it, ■ A-m ■ A)[(ij • k[) - 2 14 V] + (*4 • A)(i, ■ A-)(ij ■ *4) - 

(fcj * A*)(k 1 • A)(k2 * ^ 2 ) "1“ (^1 * ^ )(^2 * * ^ 2 ) 2^1 ^2 ] 1^1 * ^ )(^2 * * ^ 2 ) 
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+2^ V] + (k[ ■ A){k 2 ■ A*)[(fc 1 • A4) - 2fc 1 °fc'°] + {k[ ■ A)(k 2 • A*)[-( k x ■ k 2 ) + 2h% 


+(^2 • A)(k 2 • A*)p! • fc'J - 2ki°k[°] + {k’ 2 • A)(fc 2 • A^H** *1) + 2^°^°] — A - A* 

[2k x %\k 2 • fc'J - 2k 1 °k 2 °{k[ ■ k' 2 ) + k[°k 2 °{k x ■ k' 2 ) - k'^k'^ih • k&. 30) 


More generally /(ii, i 2 , *i,* 2 , *1,* 2 ) = f a ’ b (ki,k 2 ,k[,k 2 ) where A = A a andA* = A b . 


and 


g(ki,k[,k 2 , T\, t 2 , r]i,r} 2 ) = — 


(6.31) 


with 

c=4ylm° 


(6.32) 


and 

l - A - A*[2kf{k, • k[) + 2k l °k 2 °(k 2 ■ k[) - 6*1 V(fci * 2 )]+(*i ’ A)(*a • A*)(6*l V) 

+(k 2 • A)(*i • A")(6k[°k 2 0 )— (k[ ■ A){k x ■ A*){2k 2 ° k 2 °) + (k[ • A )(k x ■ A)(2k 2 °k 2 ) 

~(k[ ■ A)(k 2 ■ A’)(2ki°k 2 °) - (k[ A*)(k 2 ■ A)(2ki°k 2 °) + A°[(h • A*)2k° 2 {k 2 • k[) 

+{k 2 A*)2k 2 (ki k[)-{k[-A x )2k 0 2 (k 1 -k 2 )} + A° ,t [(k[-A)2k 0 2 {k 2 -k 1 ) 

-(h ■ A)2k 2 (k 2 • k[) - (k 2 ■ A)2kl[k 1 • *1)] (6.33) 

( t* ■ 6 1 ') 

Here A 0 = - smh^G' (t x ) , A 3 = cosh^G'in), Aq = -sinhrfrG'fa), A 3 = cosh^G' (r 2 ) and 
A 

other A’s are zero. 



Chapter 7 


J/'ty suppressions in a thermally 
equilibrating quark-gluon plasma at 
RHIC 


7.1 INTRODUCTION 

As mentioned earlier, dileptons and single photons have long been proposed as useful probes 
of the plasma[7], as once produced, they hardly interact with the strongly interacting matter 
and thus carry the details of the circumstances of their production. However, as they are also 
produced via hadronic decays m the later stages of evolution, where plasma expands and cools, 
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it becomes difficult to distinguish between them. In that sense the production of cc is a clean 
process, as it is not produced m the later stages due to its heavy mass. The production of cc 
takes place mainly at the hard vertex. Once produced, cc will evolve in to open D mesons and 
charmonium states, such as J/t/j, Xi etc., while travelling through this dense phase of matter. 
Unlike hadronic collisions, where an enhancement of J/0 is found at tevatron energy [72], both 
cc and ,//?/> get suppressed while traveling through different stages of quark-gluon plasma. 
There might be thermal charm productions in the very early stage of the plasma, where the 
temperature is very high, but no production of cc can occur in the later stages due to decrease 
in the temperature. 

We again mention that, the various stages by which the complete evolution of quark-gluon 
plasma is described in URHIC are, i) pre-equilibrium, ii) equilibrium, where one actually studies 
the thermalised quark-gluon plasma, iii) cooling and iv) hadronisation. The production of an 
equilibrated quark-gluon plasma in stage (ii) crucially depends on the pre-equilibrium evolution, 
z.e. on stage (i). Matsui and Satz[5] have suggested the suppression of J in the equilibrium 
phase. In their study, the debye screening length which is calculated from lattice QCD, is found 
to be less than the J/tp radius. This forbids the binding of cc to J/^ in the equilibrated quark- 
gluon plasma. However, such a study is not available for an equilibrating quark-gluon plasma. 
There could be many interesting effects in the pre-equilibrium stage, when cc or J/V> travel 
before reaching the equilibrium stage. As high energy deconfined partons are also present 
in the pre-equilibrium stage, J/ijj is suppressed due to its interaction with these deconfined 
partons. Since the J/U dissociation cross section inside a deconfined partons is very different 
from that inside a hadronic gas (see section II), the study of J/V> suppression can then provide 
us information on color deconfinement in the quark-gluon plasma. 
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Another outcome of the pre-equilibrium study is the prediction of the equilibration time, the time 
at which quark-gluon plasma equilibrates. A detailed knowledge of equilibration time and J/ib 
formation time is very crucial m determining the J/ip suppressions in URHIC. This is because, 
a large equilibration time may permit a J/xp formation from cc pairs before plasma equilibrates. 
Then the interaction of a fully formed J/xp with deconfined partons is prefered than screening. 

If the J /ip formation time is greater than the equilibration time or comparable to it, then cc 
interacts with the deconfined partons before screening starts operating m the equilibrium stage. 
The interaction of a cc with the partons before plasma equilibriates is not studied yet, and such 
a calculation is beyond the scope of this thesis. 

As mentioned earlier, a careful study of pre-equilibrium evolution of quark-gluon plasma is the 
consistent way to study production and equilibration of quark-gluon plasma in URHIC. From 
this point of view it is necessary to study what happens actually to the J/ip suppressions at 
different stages of evolution of QGP, rather than estimating it in an equilibrated quark-gluon 
plasma. We study here the suppression of a fully formed J/ip in equilibrating quark-gluon 
plasma using short distance QCD. The J/ip may also be suppressed in later stages due to 
inelastic collisions with hadrons, but we do not consider it here. This suppressions will be less 
due to lesser value of the hadron- J/ip inelastic cross section (described in section II). 

However, J/ip suppression is also measured m p-A collisions where there is no existence of 
quark-gluon plasma[73, 55]. At these experiments the suppression of J/ip is due to the presence 
of a nuclear medium. The prominent mechanism by which these datas are explained at SPS is 
via nuclear absorption, the suppression of a nascent J /ip before it forms a physical resonance[56, 
74]. Within this mechanism, one does not assume the existence of a quark-gluon plasma phase 
or a highly densed deconfined partonic medium. This is the case with some light nuclei collisions 
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at SPS[73, 55]. However the recent measurements of J/ip suppression by NA50 collaboration^?] 
(Pb-Pb collisions at — HGeV ) yeild an excess suppression of which is not explained by 
the above conventional approach. Presence of a nuclear medium may not be enough to explain 
this data. There are speculations that this suppression is due to the existence of a deconfined 
partomc medium[58] or a high energy density matter[59]. 

Whether an equlibrated plasma has formed in Pb-Pb collisons at SPS(v^ = 17 GeV) or not, is 
debatable, but at RHIC(Au-Au collisions at y/s = 200GeV) we may go closer to the equilibrated 
quark-gluon plasma. In this chapter we estimate the J/'tp suppressions at RHIC for a thermally 
equilibrating quark-gluon plasma. 

However before studying charmonium suppressions in quark-gluon plasma it is essential to study 
the corresponding situations in hadronic collisions. The experimental analysis of quarkonia pro- 
duction in high energy hadronic collisions at tevatron energy[72] revealed a drastic disagreement 
with the predictions from color singlet model[75]. A substantial fraction of quarkonium produc- 
tion at very high energy arises from intermediate color octet fluctuation during the evolution 
into final state color singlet Q Q pair. Now a rigorous theoretical foundation has been achived 
by Bodwin, Braaten, and Leepage[76], who devloped a QCD formalism by marrying pQCD 
and an effective field theory within non-relativistic QCD(NRQCD). This formalism accounts 
for the production of both color singlet and color octet cc states, that evolve m to final color 
singlet quarkonium. Within NRQCD, the production amplitude is expanded in powers of both 
the strong coupling a s and the velocity of heavy quark(u 2 = -23 for cc system), which includes 
higher fock-state components. The long distance non-pertubative matrix element is either fitted 
from experiments or taken, in principle, from lattice QCD[77]. This technique has been used 
to study many heavy quarkonium production processes in hadronic collisons. However a space 
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time evolution of quarkonium using this technique has not been achived. This is crucial for 
stud\ ing the evolution of J/xjj in nuclear collisions. The existence of background chromoelectric 
field [45] and high dense partonic medium complicates the study of the evolution of cc to J/^, 
both at RHIC and LHC. As cc takes 1.0 fm (time in cc rest frame) to form a J/ip bound state, 
it interacts with the nearby light partons and with the background chromoelectric field during 
this period. No attempts is made on this line to study J/'ij) production and its subsequent evol- 
ution in nucleus-nucleus collisions. This might be the consistent way to study J j^> suppression 
in URHIC. 

Since we study the suppression of J/t/i in equilibrating quark-gluon plasma, it is relevant to 
discuss the process of equilibration of QGP in URHIC, which is not completely studied yet. The 
rate of equilibration is different for diiferent models (see chapter 8). Some relevant models that 
describe the equilibration of QGP are parton cascade model(PCM)[17], heavy-ion jet interaction 
generator(HIJING)[78] and color flux-tube model. PCM and HIJING are both perturbative 
QCD based models, and color flux-tube model is a field theoritical model. The former models 
describe the evolution of hard and semihard partons and the later describes the production of 
soft partons via non-perturbative Schwinger mechanism. For a complete study of production 
and equilibration of quark-gluon plasma, color flux-tube model may be combined with the 
pQCD based models. We use here the distribution of partons from parton cascade model to 
study suppressions in a thermally equilibrating quark-gluon plasma at RHIC. A study of 
«//?/) suppressions in a chemically equilibrating quark-gluon plasma within HIJING model can be 
found in the ref. [10]. As mentioned in chapter 4 the J/t/> suppressions in a thermally equilibrating 
quark-gluon plasma within color flux-tube model will be calculated once hard partons are added 
to it via a minijet source term. Within parton cascade model, a relativistic transport equation 
is solved to study the distribution of partons by considering the direct collision processes 2 ->■ 2, 
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the inelastic collision processes 2 — >■ 1, and the decay processes 1 -» 2[17]. More explicitely, 
this is written as 

ifdufifap) = C t {x,p) l 7 - 1 ) 

where / t (x,p) is the distribution of a particular type of parton i m the usual phase space. The 
collisions term C,(x,p) is derived from pQCD taking all the processes mentioned above. The 
initial condition on f t (x,p ) is the distribution of partons inside the nucleus before the two nuclei 
collide with each other [38]. However for a complete study, the addition of soft partons and the 
effect of background field has to be taken into account. These backgruond electric field may 
not exit in the very early stage where pQCD is applicable[41], but in later times this will have 
an important role in the equilibration of the plasma[41, 40]. Such calculations may be available 
once the color flux-tube model is combined to the pQCD based models. 

In section II we describe the short-distance QCD and the dissociation of Results and 

discussions are presented in section III, for RHIC energy. We summarise and conclude the 

main results in section IV 


7.2 J / ip dissociation by deconfined gluons 

In the framework developed by Bhanot and Peskin[79, 80], interactions between light hadrons 
and deeply bound, heavy quarkonmm states, such as the -W, axe mediated by short-range 
color dipole interactions. For sufficiently heavy quarks, the dissociation of quarkonium states 
by interaction with light hadrons is fully accounted for by short-distance QCD[79, 80]. Because 
of its small size, a heavy quarkonium can probe the short distance properties of light hadron. A 
parton based calculation of the W-hadron cross section is thus possible via an operator product 
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expansion method, simillar to that used in deeply inelastic lepton- hadron scatterings [79-82]. 
These pertubative calculations become valid when the space and time scale associated with 
the quarkonium state, tq and tq, are small in comparison to the nonperturbative scale 
which is the characteristic size of the light hadrons, i.e. tq «< A an< ^ * Q <<< ^qcd- 
For J/i/> ground state rj ^ ~ 0.2/m = (1 GeV)~ l and Ej/^(2Md — M^) ~ 0.64 GeV. With 
A qcd — 0.2 GeV , the above inequalities seem to be well satisfied[54] and one expect that the 
dissociation of <//?/> by hadron will be governed by the J/^- hadron break-up cross section as 
calculated m short distance QCD. The operator product expansion allows one to express the 


hadron -J/ip inelastic cross section in terms of the convolution of the gluon- J/?/ dissociation 
cross section with the gluon distribution inside the hadron. The gluon- J/ip dissociation cross 
section is given by[79, 80, 82]: 


2rr 


. 1 6 7T , 


( O', _ ~ n (OC) VI !m 2 ', ( < ?°/ £ ° ~ 1 ) 3/ ~ 

~ J^ 32 / 3) ( 3g2 (q°/e o y ' 


(7.2) 


Here g 3 is the coupling between gluon and charm quark, m Q the charm quark mass, and q° the 
gluon energy in the Jji> rest frame. The heavy quarkonium is coulomb like, and one writes, 


eo = (f&W 79, 82], 


The gluons which are soft inside a pion are not capable of dissociating a charmonium. Hence 
they give a low ■n-J/tp cross-section. However, deconfined gluons, which carry enough energies, 
are sufficient to break a charmonium [54]. These high energy partons are present in DRHIC, 
as observed in several calculations, such as[45, 17]. These deconfined gluons can break a fully 
formed J/ip that exit inside QGP. This conclusion does not seem to be affected substantially 

by nonperturbative effects [83]. 

We use the above g-J/i/, dissociation cross section to study the J hb suppressions in a thermally 
equilibrating quark-gluon plasma at RHIC. For the central collisions, the J/i, survival probab- 
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llity is: 


S 9-J/ip( P t) = exp[- J drn g (r) < v rel cr(k • v) >*]. (7.3) 

Heie n g [r) is the gluon number density which evolves as system expands longitudinally. Accord- 
ing to Bjorken sceneno[2], this number density is a function of the boost invariant parameter r 
( r ~~ V^ 2 ~~ thermal average g -J/4> cross section, < v re i<j(k • v) >*., is 




(7.4) 


Here v(= (Mt,Pt,0)/ Mj/^) is the four-velocity of Jjip in central rapidity region and k is the 
four-momentum of gluon in the parton gas. In the thermal average cross-section we use the dis- 
tribution function, /(fc°,T(r)) = e X p( A a? ( T )) — 1 ^ or S^ uon - Here a(r ) captures the deviation from 
equilibrium, which m principle, is determined from the relation, n g (r ) = f dT(p tl u^)f(p,T(T)), 
by knowing the evolution of temperature and n g , for an equilibrating quark-gluon plasma. Here 
dr = - ^ 3 ^- with 7 = 2 x 8 (the product of spin and color degeneracy), and u( = (t/r, 0 , 0 , 
z/r)) is the flow velocity. However, as a(r) cancels both from numerator and denominator 
of equation(7.4), we need not know its form. What matters is n 3 (r) and T(r) whose evolu- 
tions are taken from PCM[17]. Within PCM, n g (r) ~ 0.72n(r), with n(r) = 565/m~ 3 (^-)~°' 90 


and T(r) = 950iVfeV( 


X '1-0.30 

V 5 


for a thermally equilibrating quark- gluon plasma at RHIC. Here 


To = .05 fm. Using this distributions the integration in equation (7.3) is performed numerically 
to study the J /ip survival probability. 


7.3 Results and discussions 

In Fig- 7.1 we present the thermal averaged gluon -J/ip cross section as a function of J/rp trans- 
verse momentum Pj, for different values of temperature. It can be seen that this value is de- 
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creased as the temperature of the plasma is increased. Also for J/ip with high P T this thermal 
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Figure 7.1: The thermal averaged gluon- J/tp cross section < v re \a > (mb) as a function of 
transverse momentum Pj at different temperatures. Solid line refers to T= 0.2 GeV, upper 
broken line corresponds to T= 0.4 GeV and lower broken line refers to T=0.8 GeV. 

cross section becomes less. This has a direct impact on J /ip suppression. As Pt becomes large 
the survival probability must increase. However, due to large initial parton density formed m 
the initial stage of the plasma, these smaller thermal cross sections do not have a large effects 
on the total survival probability. This is shown in Fig-7.2 where we have plotted the survival 
probability S g -j/^ as a function of Pt for a thermally equilibrating QGP at RHIC within PCM. 
As can be seen, the suppression is almost 100 percent for J/ip with low transverse momentum. 
Even for J/ip with high Pt(— 10 GeV) the suppression is about 95 percent. This huge suppres- 
sion is due to the presence of a highly densed deconfined partonic medium at RHIC. As can 
be seen from equation(7.3) and Fig-7.1 this suppression will be more for a plasma with larger 
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Figure 7.2: The survival probability of J/ip in a thermally equilibrating plasma at RHIC. 


number density and lower temperature. A denser and cooler plasma is obtained once the soft 


partons studied within the color flux-tube model in earlier chapters are added to pQCD based 


models such as PCM. The measurement of such huge suppression will reveal the existence of 


such a deconfined partonic medium, possibly a thermalised quark-gluon plasma, at RHIC. 


However, what one actually measures experimentally is the total survival probability This 


implies that other sources of J/ib suppression, without the existence of a QGP phase, have to 


be singled out from the total survival probability in order to see the suppressions only from this 


deconfined partonic medium. If the estimated suppression due to all sources without QGP phase 


is much less than the experimental findings, one may then hope that the rest of the suppressions 


is due to the existence of a highly densed deconfined partons. The experimental measurement 


of J jib at RHIC will be able to reveal all the possibilities and may explain the existence of such 


a deconfined partonic medium, possibly a thermalised quark-gluon plasma. 


On the other hand J/i/> enhancement can occur due to an increase in the (thermal) charms 





99 


produced at a very high temperature, especially in the initial stage of the evolution of the 
system [bo, 36] We recall that screening of J ftp was origmally[5] proposed without taking 
production of thermal charms into account, which was argued to be negligible at T ~ 200 - 300 
iVIeV However, the system may approach a temperature which is as large as 900 MeV initially, 
and can therefore produce a large number of thermal charms at RHIC So to make resonable 
estimates at RHIC and LHC, one has to take into account both thermal and hard J/tp. For 
simplicity, we have not considered here those J/ip which are created from thermal charms. 


7.4 Conclusion 

Applying short distance QCD mechanism and taking parton cascade model for the evolution of 
quark-gluon plasma at RHIC we have found a huge suppression (almost 100 percent) of J/ip 
which are produced from primary charm quarks. In principle, the total survival probability is 
written as 

St — Sg^J/^Sother ( 7 - 5 ) 

where first one ( S g -j /$) is due to g -J/ip dissociation, and the last on e{S ot her) is due to any other 
possible sources of suppression^ 1-14] such as nuclear absorption, final state comover scatterings 
etc. As the hadron- J/xf) cross section is smaller than g -J /ip cross section[54], we expect a lesser 
suppression of J/ip in the hadromsation stage. The prominent suppression seems to be due to 
the nuclear absorption (60 percent at RHIC[9]) and J / 'p dissociation m the deconfined partonic 
system. Hence the measurement of a huge Jjip suppression at RHIC may suggest the existence 
of a partonic plasma. However, many refinements, such as space time evolution of quarkonium 
production m nucleus-nucleus collision using NRQCD has to be worked out, before unambigous 
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conclusions can be drawn for the pre-resonance suppression of J ftp. Also, the contribution from 
thermal .7/0, which may not be neglible at RHIC and LHC, has to be taken into account. It is 
only after incorporating the above features the certainty of screening can be justified. 

In any case, we have shown here the importance of the thermally equlibrating quark-gluon 
plasma on J/0 suppression. This suppression depends on variuos models which describes the 
thermal evolution. From the measurement of this suppression, the validity of these models for 
thermal evolution of quark-gluon plasma at RHIC can be checked. 



Chapter 8 


A brief description of other Models and 
approaches 


8.1 Dynamical Models in heavy-ion collisions 

8.1.1 Hevy-ion jet interaction generator (HIJING) 

H eaV y_ion jet interaction generator model[37] is a pQCD inspired model where the non-perturbative 
part is modeled phenomenologically The evolution of hard and semihard partons are described 
by pQCD. The mtitial condition is estimated by HIJING monte carlo code in which initial dir- 
ect parton scatterings are taken into account. The produced parton gas in the central rapidity 
region approaches thermal and chemical equilibrium through elastic scatterings and induced ra- 
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diations. In an overlap region of lfm thickness (of the two nuclei) a local isotropy in momentum 
distributions assumed and is estimated to be achieved in a time period of 0 . 5 - 0 . 7 fm through 
free-streaming evolution. It is further assumed that this is the actual kinetic equilibration ( or 
thermal isation) time for the partonic system. Though this assumption is not well founded, this 
is the best that, the model incorporates in the absence of a complete calculation of the relaxa- 
tion time. Within this model the scattering processes which maintain thermal equilibrium are 
,jg 44 gg, g q <4 gq etc.. Further chemical equilibration of partons is studied by by solving master 
rate equations, where gg ggg and gg qq processes are taken into account. Asssummg that 
the parton scatterings are sufficiently rapid so as to mainatin thermal equilibrium, and therefore 
neglecting the effects of viscosity due to elastic and inelastic scatterings, one determines the 
evolution of the energy density from the hydrodynamic equation. The master equations are 
solved self consistently[37] to obtain the time evolution of the temperature and fugacities. It is 
found that the parton gas cools considerably faster than the Bjorkens scaling law (T 3 r = cons.), 


because the production of additional partons during the approach to chemical equilibrium con- 
sumes an appreciable amount of energy. It is shown within this model that the equilibration 
can hardly reach before the temperature drops below 200 MeV in the short period time of 1-2 
fm at RHIC energy. However, at LHC energy the parton gas, especially the gluon component, 
becomes very close to equilibrium and the plaama may exist in a deconfined phase for as long 
as 4-5 fm. One observation of this model is that the quarks never reach chemical equilibrium at 


both energies. 
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8-1.2 Parton Cascade Model (PCM) 

I he parton cascade model[17] is a relativistic kitetic model for high energy nuclear collisions 
which is inspired by the parton picture of hadronic collisions. The dynamics of the space time 
evolution of quarks and gluons is studied by the interaction of partons m the framework of 
perturbative QCD. The time devlopment of parton cascades is taken in six-dimensional phase 
space of coordinate and momenta according to a semiclassical, relativistic transport equation 
with a collision kernel that contains the asymptotically dominant perturbative QCD interaction 
processes among the partons. The time devlopement of the parton distributions, the n um ber 
and type of interactions is determined by the dynamics itself. The transport equation which is 
solved is 

/a,F.(r, f ) = SJimr.f) (8.1) 

where F a (r,p) is the phase-space distribution of the parton of type a and l[k) is the collision 
kernel for the collision terms. The initial state for the parton phase space distribution is taken 
from the experimantally measured nucleon structure functions and elastic nucleon form-factors, 
as well as nuclear shadowing corrections. The general form of collision term, cr k l[K)(r,p), in the 
above equation includes a subset of all possible n -> n' interaction proceses. These processes 
consists of 2 -+ 2 colisions, 2 -> 1 fusions and 1 -» 2 decays of partons in the lowest order 
perturbative QCD, plus higher order corrections which are associated with multiple emissions 
and absorptions of partons. 

Within this study it observed that isotropy in the momentum distribution is reached very 
quickly, leading to very rapid thermalisation, almost after 0.2 fm of maximum overlap of two 
nuclei at RHIC. The system is very dense (n(r) = 565 /m- 3 ^)' 090 ) and very hot (T(r) = 
%0MeV(^)~° 90 ) within this study. 
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8.2 Lattice QCD studies 

In the high temperature limit of QCD, the short range interactions among the partons become 
weak and the long range interactions become dynamically screened. As a consequence, hadronic 
matter at very high temperature do not exhibit confinement and chiral symmetry breaking. 
Lattice gauge theory which allows nonperturbative numerical calculation of observables m QCD 
predicts that the transition between low-temperature and high-temperature phases of QCD is 
not smooth but discontnuous. The order of chiral phase transition is sensitive to the number 
of light, dynamical quark flavors. This phase transition is of second-order for two massless 
flavors and is of first-order for three massless flavors. The numerical simulations of lattice gauge 
theory have predicted that the transition temperature lies m the range 150±10 MeV at vanishing 
net quark density. The simulations without dynamical quarks clearly exhibit a first-order phase 
transition at a temperature T c ~ 260 MeV. The quarks are not confined in the high temperature 
phase. The quark condensate < ^ > shows a strong drop over the same temperature range, 
indicating that deconfinement and restoration of chiral symmetry go hand-in-hand. However, 
the simulations of lattice QCD with dynamical quarks have not yet overcome the limitations due 
to finite lattice size. More works in this area are now in progress regarding the order of phase 

transition. 



Bibliography 


[1] J. W. Harris and B. Muller, in Annual Review of nuclear and particle sciences, Vol. 46 
71 (199C) edited by Chris Quigg al. . 

[2] .J. D. Bjorken, Phys. Rev D27 (1983)140 

[3] S. E. Eiseman et a/., Phys. Lett. B297, 44 (1992). 

J. Barrette et al., Z. Phys. C59, 211 (1993). 

T. Abbott et al., Phys. Rev. C50, 1024 (1994). 

S. E. Eiseman et al., Phys. Lett B325, 322 (1994). 

J. Barrette et al., Phys. Lett. B351, 93 (1995). 

[4] E. Andersen et al., Phys. Lett. B327, 433 (1994). 

M. Gadzicki et al., Nucl. Phys. Lett. A590, 197c (1995). 

D. D. Bari et al., Nucl. Phys. Lett. A590, 307c (1995) and refernces therein. 

[5] T. Matsui and H. Satz, Phys. Lett . B178 416 (1986) 

[6] M. T Strickland, Phys. Lett. B331 245 (1994). 

[7] Jan-e Alam, Sibaji Raha and Bikash Sinha, Phys. Rep. 273 243 (1996). 


106 



107 


[8] J. Rafelski and B. Muller, Phys . Rev. Lett. 48 1066 (1982), C. P. Singh, phys. Rev. 
Lett. 56 870 (1990). 

[9] Gouranga C Nayak, Journal of High Energy Physics 02 (1998) 005. 

[10] Xiao-Ming Xu e.t a/., Phys. Rev. C 53, 3051 (1996). 

[11] S. Gavin and R. Vogt, Phys. Rev. Lett. 78, 1006 (1997). 

[12] C. Y. Wong, Phys. Rev. C 55, 2621 (1997). 

[13] R C Hwa, T. Pisut and N. Pisutova, OITS 628 (1997). 

[14] Rajiv V Gavai and Sourendu Gupta, TIFR/TH/97-07, Report No:- hep-ph/9703300, To 
appear m Phys. Lett. B. 

[15] J. Rafelski and R. Hagedorn, in statistical Mechanics of Quarks and Hadrons, Ed. H. 
Satz, (North Holland, Amsterdam, 1981). 

[16] S. Nagamiya, Nucl. Phys. A544 5c (1992). 

[17] K. Geiger and J. Kapusta, Phys Rev D47 4905 (1993), and references therein. 

[18] G. Baym, Phys. Lett. B138 18 (1984). 

[19] K. Kajantie and T. Matsui, Phys. Lett B164 373 (1985). 

[20] B. Andersson et al. Phys. Rep. 97, 31 (1983), Nucl. Phys B 281, 289 (1987). 

[21] F E. Low, Phys. Rev. D12 163 (1975). 

[22] S. Nussmov, Phys. Rev. Lett. 34 1286 (1975). 

[23] J. Schwinger, Phys. Rev. 82 664 (1951). 



108 


[24] A Bialas et el. Nucl. Phys. B296 611 (1988). 

[25] B. Banerjee, R. S. Bhalerao and V Ravishankar, Phys. Lett. B224 16 (1989). 

[26] T T. Wu and C. N. Yang, Phys. Rev. D12 3843 (1975). 

[27] R V Gavai, private communication. 

[28] M. Asakawa and T. Matsui, Phys. Rev. D43 2871 (1991). 

[29] (J. Heinz, Phys. Rev. Lett. 51 351 (1983). 

[30] S. K. Wong, Nuovo Cemento A65 689 (1970). 

[31] H-T Elze and U. Heinz, Phys. Rep. 183 81 (1989) 

[32] U. Heinz, Ann. of Phys. (NY) 161 48 (1985). 

[33] R. Pisarski, Phys Rev. Lett. 63 1129 (1989). 

[34] P. F. Kelly, Q. Liu, C. Lucchesi and C. Manuel, Phys. Rev Lett. 72 3461 (1994) 

[35] P. F. Kelly, Q. Liu, C. Lucchesi and C. Manuel, Phys. Rev. 50 4209 (1994) 

[36] B. Svetitsky, Phys. Rev D37 2484 (1988). 

[37] X. N. Wang, Phys. Rep. 280, 287 (1997). 

[38] K. Geiger, Phys. Rep. 258, 237 (1995). 

[39] R. S. Bhalerao and V. Ravishankar, Phys. Lett B 409, 38 (1997). 

[40] K. .J Eskola, B. Muller and X. N. Wang, Phys. Lett. B 374, 20 (1996). 

[41] K. J. Eskola and M. Gyulassy, Phys. Rev. C 47, 2329 (1993). 



109 

[42] L. S. Brown and W. I. Weisberger, Nucl. Phys B. 157, 285 (1979). 

[43] H. Heiselberg and X. N. Wang, Phys. Rev. C 53, 1892 (1996). 

[44] C M. Hung and E.V Shuryak, report no hep-ph/9709264 

[45] Gouranga C. Nayak and V. Ravishankar, Phys. Rev. D55, 6877 (1997): and references 
therein. 

[46] Gouranga C Nayak and V. Ravishankar, Phys. Rev C58 356 (1998). 

[47] M. Gyulassy and A. Iwazaki, Phys. Lett. B 165B, 157 (1985). 

[48] M. Gyulassy and L. P. Csernai, Nucl. Phys. A 460 723 (1986). 

[49] Increasing the accuracy by another 5 percent does not substantially alter the estimate. 

[50] G. C. Nayak and V. Ravishankar, in preparation. 

[51] see ref[41]. The Ohmic heating in the reference is realised by the particle production in 
our model. 

[52] Stanislaw Mrowczynski, Phys. Lett. B 393 26 (1997). 

[53] see, e.g., K. Rajagopal, CALT-68-2104, who has discussed the possible impact of non- 
equilibrium state on the formation of disoriented chiral condensates. 

[54] D. Kharzeev and H. Satz, in Quark-Gluon plasma II, edited by R. C Hwa (World Scientific 
, Singapore, 1995), p. 395 

[55] C. Baglin et al. Phys Lett. B 270, 105 (1991); 


C. Baglin et al. ibid 345, 617 (1995). 



110 

[56] C. Gerschel and J. Hufner, Z. Phys. C 56, 171 (1992). 

[57] M. Gonin et al. Nucl Phys. A 610, 404c (1996). 

[58] D Kharzeev, C. Lourenco, M. Nardi and H. Satz, Z. Phys C 74, 307 (1997) 

[59] J. P. Blaziot and J. Y. Olhtrarault, Phys. Rev. Lett. 77, 1703 (1996) 

[60] R. C Hwa, T. Pisut and N. Pisutova, OITS 628 (1997) 

[61] B. Kampfer and O. P. Pavlenko, Phys. Lett. B289 127 (1992). 

[62] L. D. McLerran and T. Toimela, Phys. Rev D 23 545(1985) 

[63] K. Kajantie, M. Kajata, L. McLerran and P. V. Ruuskanen, Phys Rev. D 34 811 (1986). 

[64] A. Bialas and J. P. Blaizot, Phys. Rev. D 32 2954 (1985). 

[65] K. Geiger and J. I. Kapusta, Phys. Rev. Lett 70 1920 (1993). 

[66] J. Kapusta, L. McLerran and D. K. Srivastava, Phys. Lett. B283 145 (1992). 

[67] A. Karman, T. Matsui and B. Svetitsky, Phys. Rev. Lett. 56 219 (1986) and references 
therein. 

[68] A. Casher, H. Xeuberger and S. Nussinov, Phys. Rev. D20 179 (1979). 

[69] Coherent states, Applications in Physics and Mathematical Physics, edited by J. R. 
Klauder and Bo-Sture Skagerstam (World scientific, 1985). 

[70] Mark Byrd, DOE-ER-40757-104, UTEXAS-HEP-97-18 (1997), Report No:- phys- 
ics/9708015. 



Ill 


(/ 1] A. K. Ganguli, P. K. Kaw and J. C. Parikh, Phys. Rev. D48 R2983 (1993) and references 
therein. 

[i'2] M. Mangano, proceedings of the X’th Topical Workshop on pp physics, Fermilab public- 
ation 1995(hep-ph/9507353). 

[73] ,J. Badier et al. Z Phys. 20, 101 (1983); D. M. Aide et al. Phys. Rev. Lett. 66, 133(1991). 

[74] D. Kharezeev and H. Satz, Phys. Lett. B 366, 316 (1996). 

[75] R. Baier and R. Ruckl, Z. Phys. C 19, 251 (1992). 

[76] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D. 51, 1125 (1995). 

[77] Peter. Cho, and Adam. K. leibovich, Phys. Rev. D 53, 150 (1996). 

[78] T. S. Biro, E. Van Doom, B. Muller and X. N. Wang, Phys. Rev. C 48, 1275 (1993). 

[79] M. E. Peskin, Nucl. Phys. B 156, 365 (1979). 

[80] M. E. Peskin and G. Bhanot, Nucl. Phys. B 156, 391 (1979). 

[81] A. Kaidalov, in QCD and High Energy Hadronic interactions, edited by J.Tran thanh 
Van (Edition Frontier, Gif-sur- Yvette, 1993) 

[82] D. Kharezeev and H. Satz, Phys. Lett. B 334, 155 (1994). 

[83] D Kharezeev, L. McLerran, and H. Satz, Phys. Lett. B 356, 349 (1995). 

[84] C Y. Wong, Phys. Rev. Lett. 76, 196 (1997); C. W Wong, Phys. Rev. D 54, R4199 
(1996). 


[85] K. Geiger, Phys Rev. D 48, 4129 (1993). 



112 


[86] A. Shor, Phys. Lett. B 246, 207 (1990). 



